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Abstract 

Quantum transients are temporary features of matter waves before they reach a stationary 
regime. Transients may arise after the preparation of an unstable initial state or due to a sudden 
interaction or a change in the boundary conditions. Examples are diffraction in time, buildup 
processes, decay, trapping, forerunners or pulse formation, as well as other phenomena recently 
discovered, such as the simultaneous arrival of a wave peak at arbitrarily distant observers. The 
interest in these transients is nowadays enhanced by new technological possibilities to control, 
manipulate and measure matter waves. 
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1. Introduction 



Transients are temporary phenomena before a system attains a steady-state responding 
to a sudden change, such as an interaction switch, an excitation, new boundary condi- 
tions, or the preparation of an unstable initial state. They are ubiquitous in physics, and 
frequently described by characteristic times: response times, build-up times, decay times, 
or times of emission and arrival. In quantum mechanics, matter-wave transients are pe- 
culiar and generally distinct from corresponding classical particle phenomena. Typical 
transients are damped oscillations, forerunners, response to excitations, and trapping or 
decay processes. 
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In this work we shall review the theory, results, and experiments of quantum transients 
of single to few-body systems. We shall leave aside many-body and complex systems, such 
as electrons in a semiconductor responding to a strong short optical pulse, since they are 
treated already in several monographs (see e.g. [l|2j) and require specific techniques. 
We shall nevertheless discuss many-body systems in which effective one-particle pictures 
are still a useful approximation, e.g. a Bose-Einstein condensate within the mean field 
Gross-Pitaevskii equation, Tonks- Girardeau gases, or multi-electron systems treated at a 
Hartree level. Many of the concepts and techniques discussed here are valid for studying 
exponential-decay processes and deviations from it, but we shall not include decay, or 
touch it only in passing, since the amount of work on that field and its applications deserve 
separate monographs and many aspects have already been reviewed "Sl. The field covered, 
in spite of the above exclusions, is still vast, and we have paid special attention to work 
done by the authors, while making an effort to offer a global perspective. 

Quantum transients have frequently been studied per se, because of their interesting 
dynamics. They are often analytically solvable and this has also motivated their use to 
clarify fundamental questions such as tunneling dynamics, or the Zeno effect. In recent 
times, with the advent of technological advances such as laser cooling techniques, meso- 
scopic devices, ultrashort laser pulses and nanotechnology, experiments and applications 
which were unthinkable one or two decades ago can be implemented or are within reach. 
This is the case of operations to prepare and control the dynamics of quantum systems for 
fundamental studies, interferometry and metrology, or for quantum information process- 
ing. An understanding of transients is therefore essential: sometimes just to avoid them, 
as in an atom-laser beam; but also to take advantage of them, for example to transmit 
information simultaneously to receivers at unknown distances; to optimize gate opera- 
tions or transport properties in semiconductors, or as diagnostic tools for determining 
momentum distributions, interactions, and other properties of the system. 

A common thread in the mathematical description of many quantum transients is the 
Moshinsky function, which is related to the Faddeyeva w function, to Fresnel integrals 
and to the error function. Much of this review is devoted to study phenomena where 
the Moshinsky function plays a fundamental role. It was introduced with the method 
of time-dependent boundary conditions used in the dynamical description of heavy-ion 
collisions [4]. It was later associated with the "quantum shutter" problem, one of the 
archetypal quantum transient phenomena, which lead to the concept of "diffraction in 
time" (DIT) [5]. A flurry of experimental and theoretical work has been carried out ever 
since. The hallmark of DIT consists in temporal and spatial oscillations of Schrodinger 
matter waves released in one or several pulses from a preparation region in which the 
wave is initially confined. The original setting consisted in a sudden opening of a shutter 
to release a semi-infinite beam, and provided a quantum, temporal analogue of spatial 
Fresnel-diffraction by a sharp edge [5j. Later on, more complicated shutter windows 
[6], initial confinements |7l8l9ll0j . time-slit combinations [llll2j and periodic gratings 
[13|14|15] have been considered [16]. The effects of an external linear potential, e.g. due 
to gravity, on the DIT phenomenon have also be taken into account [17] . The basic results 
on Moshinsky 's shutter are reviewed in Section [2l 

The quantum shutter problem has important applications indeed, for example the mod- 
elization of turning on and off a beam of atoms, as done e.g. in integrated atom-optical 
circuits or a planar atom waveguide [18] , and has been used to translate the principles of 
spatial diffractive light optics |19|20j to the time domain for matter waves [2Tj . It is at the 
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core of different schemes for loading ultracold atoms or Bose-Einstein condensates into 
traps [2? and, besides, it provides a time-energy uncertainty relation [23124125] which 
has been verified experimentally for atomic waves in a Young interferometer with tem- 
poral slits |26|27|28j . The Moshinsky shutter has been discussed in phase space [29|30) . 
for relativistic equations |5|31|32|33|34|35|36|37j . with dissipation [38|36j . in relation to 
Feynman paths [39] , or for time dependent barriers [6ll6|40l41|42j . 

Transients for the expansion dynamics from box-like traps have also been examined 
both at single [814319] and many-particle level 4 4145146147148] . as well as from spherical 
traps similar to actual experimental magneto-optical traps [10) . The experimental buildup 
and depth control of box-like (square well) traps for ultracold atoms with an all-optical 
implementation [49] , or in microelectronic chips [5(T has excited a great deal of attention 
for these geometries. The transient effects when the trap is turned off provide information 
on the initial state and allow for its reconstruction [51 . Once the ballistic dynamics 
dominates the expansion, the regime of interactions or the momentum distribution can 
be inferred in a Bose-Einstein condensate [52|53j and effectively one-dimensional gases 
as well [54|55] . The involved time scales play therefore a key role and are reviewed, at a 
single-particle level, in Section [31 and with particle interactions in Section O 

Section 131 deals with pulse formation, a key process in ultracold neutron interfcrometry, 
and a timely subject due to the possibilities to control the aperture function of shutters 
in atom optics, the ionization by ultrashort laser pulses, or the development of atom 
lasers. Some of the first mechanisms proposed explicitly for building atom lasers implied 
periodically switching off and on the cavity mirrors, that is, the confining potential of 
the lasing mode. As an outcome, a pulsed atom laser is obtained. Much effort has been 
devoted to designing a continuous atom laser, whose principle has been demonstrated 
using Raman transitions [56|57j . With this output coupling mechanism, an atom initially 
trapped undergoes a transition to a non-trapped state, receiving a momentum kick due 
to the photon emission. These transitions can be mapped to the pulse formation, so 
that the "continuous" nature of the laser arises as a consequence of the overlap of such 
pulses. The fields of coherent control, femto and attosecond laser pulses and ultracold 
matter are also merging rapidly creating new opportunities for inducing, studying, and 
applying quantum interference of matter pulses and the corresponding transients. There 
is, in summary, a strong motivation for a thorough understanding of matter-wave pulse 
creation, even at an elementary single-particle level. 

The shutter model, with adequate interaction potentials added, see Section (5] has 
been used to study and characterize transient dynamics of tunneling matter waves 
[58|171 59 60 61 62 63 ,32] , reviewed in Section |7l Section |8| is devoted to the transient 
response to abrupt changes of the interaction potential in semiconductor structures and 
quantum dots J64'65'66]. The study of such transients has been further proposed for 
potential barrier classification [67] . 

In Section [9| some transient phenomena different from DIT which have been recently 
discovered or studied are reviewed. They are intriguing effects such as: ultrafast prop- 
agation of a peak which arrives simultaneously at arbitrarily distant observers in an 
absorbing medium; breakdown of classical conservation of energy in quantum collisions 
and momentum-space interferences; transients associated with classically forbidden mo- 
menta, e.g. negative incident momenta that affect the transmitted part of a wavepacket 
impinging on a barrier from the left; and the Zeno effect in a time-of- arrival measurement 
model. 
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The experimental work is reviewed in Section [TU) The first field where transient effects 
and a corresponding time-energy uncertainty relation were investigated was ultracold 
neutron interferometry [68169170171121172] . Experimental observations of transients have 
been carried out with ultracold atomic waves where diffraction, phase modulation and 
interference were studied using temporal slits [26|27|28] . More recently, a double tem- 
poral slit experiment with electrons has been performed |73j . The applications of atom 
interferometry in time domain have also been explored with Bose-Einstein condensates 
at both theoretical [75 and experimental level [75]. A very recent and relevant result 
is the realization of an electromagnetic analogy allowing to study DIT phenomena by 
means of the vertical propagation of the lasing modes emitted from the laterally confined 
cavities [76]. Note also that transients do not necessarily involve spatial translation and 
release from traps or shutters. They may be observed in the evolution of internal state 
populations of an atom due to time-dependent laser excitations [77] . 

The review ends up with some conclusions and technical appendices. Most sections 
could be read independently although there are strong links between Sections |3] and \E\ 
or between Sections |6] and [T] Multiple use of some symbols has been unavoidable and the 
notation is only guaranteed to be consistent within a given section. 

2. Diffraction in time 

2.1. The Moshinsky shutter 

In 1952, Moshinsky tackled the quantum dynamics of a suddenly released beam of 
independent particles of mass m. The corresponding classical problem, that of a beam 
with a well defined incident velocity p/m, and initially confined to the negative semi- 
axis (a; < 0), admits a trivial solution: the density profile is uniform for x < pt/m 
and vanishes elsewhere. In non-relativistic quantum mechanics, we shall first consider a 
quasi- monochromatic atomic beam of nominal momentum hk (or velocity Vk — hk/m) 
impinging on a totally absorbing shutter located at the origin a; = 0, 



where 9(x) is the Heaviside step function defined as Q{x) — 1 for a; > and elsewhere. 
Note that, in spite of the prominent role of fc, the state is not really monochromatic 
because of the spatial truncation. Such state is clearly not normalized, but it can be 
considered as an elementary component of a wave packet truncated (fully absorbed) at 
the origin. Since for i > the shutter has been removed, the dynamics is free and the 
time evolution can be obtained using the superposition principle, 



^(a;,t = 0) = e^'^^e(-a;). 



(1) 



OO 




(2) 



— OO 



with the free propagator 




Fig. 1. Snapshot of the density profile of a suddenly released beam, confined transversely to a harmonic 
trap along a waveguide, and exhibiting difi^raction in time in the longitudinal direction, see Eq. 

see Appendix [X] for a brief account of propagators and their role in quantum transients. 
Following customary practice, we shall liberally speak of \tp{x,t)\'^ as a "density" even if 
it is dimensionless. 

In a seminal paper [5] , Moshinsky proved that 

il^{x,t) = M{x,k,t). (4) 

This is a simple compact result, which constitutes the main reference case for more 
complex setups. M is the so-called Moshinsky function |4|5|24j . 

, r/ , ,\ , , l + i [ht f , mx\ ... 

Mix,k,t):=^M-u), -=^i-[k-^)- (5) 

Many of its properties relevant to the study of quantum transients are listed in Ap- 
pendix [BJ Here we just note that it can be related to the Faddeyeva function w{z) := 
erfc(— iz) [781 79j . see Appendix [Cl Such solution entails the diffraction in time phe- 
nomenon, a set of oscillations in the beam profile (see Fig. [T]) in dramatic contrast with 
a classical monochromatic and homogeneous beam density Q(pt/m — x). 

The "diffraction in time" (DIT) term was introduced because the quantum density 
profile, \M{x,k,t)\'^, admits a simple geometric interpretation in terms of the Cornu 
spiral or clothoid, which is the curve that results from a parametric representation of the 
Fresnel integrals C and S, see Fig. [2] and Eq. (jC.6p . As shown also by Moshinsky [S], it 
follows that 



with 
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8(9) 




C(9) 



Fig. 2. Cornu spiral. The density profile characteristic of the diffraction in time, Eq. is given by half 
the distance from the point (—1/2, —1/2) to any point of the spiral. 

This is precisely the form of the intensity profile of a light beam diffracted by a semi- 
infinite plane [19j . which prompted the choice of the DIT term. The probability density 
is half the distance from the point (—1/2, —1/2) to any other point of the Cornu spiral. 
In such representation the origin corresponds to the classical particle with momentum 
p released at time t = from the shutter position. Thanks to it, the width of the main 
fringe, Sx, can be estimated from the intersection between the classical and quantum 
probability densities [5|12j , leading to a dependence of the form 



6x cx (nht/m)^^^. 



(8) 



Moreover, any point of constant probability deviates from the classical spacetime path. 
In particular, the principal maxima and minima obey 

nht „ 



(9) 



Xmax (t) = pt/m - \ 



Xmin(t) — pt/m \ Qmin: 

m 



where 9max = 1-217, 9min = 1.872 are specific values for the totally absorbing shutter. 
Therefore, the quantum signal is slowed down for t < with respect 

to the classical case, Xdit) = pt/m. 

A more general type of initial state was considered in [24112] , 



^{x,t = 0) = e''''' + Re' 



(10) 



with R = exp(zQ!7r) corresponding to a shutter with reflectivity = 1. Under free 
evolution, tp{x,t) = M{x,k,t) + RM{x,—k,t). For such a beam, the visibility of the 
fringes increases monotonically from a = (cosine initial condition, see 17.21 below) to 
a — 1 for which the shutter is equivalent to a totally-reflecting hard-wall potential, see 
Fig. [3] and Appendix [Dl 
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Fig. 3. Enhancement effect of the reflectivity of the shutter on the visibiUty of the diffraction-in-time 
fringes of a suddenly released quasi- monochromatic matter- wave beam (in h = m = I units). 

The possibility of closing the shutter after a time t to form a pulse was studied by 
Moshinsky himself [23l24j . This is the essential step of a chopper in an ultracold neutron 
interferometer used for velocity selection. For small values of r a time-energy uncertainty 
relation comes into play, broadening the energy distribution of the resulting pulse. For 
the initial condition (fTO|) . with R = +1, the explicit calculation of the energy spread AE 
in the resulting wavepacket leads to 



Using the full width at half-maximum (FWHM) of the energy distribution, it was found 
that 



The time-energy uncertainty relation lacks the simple status of that holding for position 
and momentum observables, having many different non-equivalent versions reviewed in 
[25] . The one for chopped beams as discussed here [23|24j was experimentally observed 
and will be further discussed in Section [TOl 

2.2. Time- dependent shutter 

The paradigmatic result by Moshinsky rests on the sudden approximation for the 
removal of the shutter. The effect of a more realistic time dependence has also been 
looked into. In [6|16j . the shutter potential V{x, t) = ^5{x)/t was assumed, which reduces 
to an infinite wall at i = and vanishes asymptotically for t — > oo. Using the Laplace 
transform method, the associated propagator can be related to the free one [6j, from which 
the time-dependent wavefunction can be calculated by numerical integration. Scheitler 
and Kleber analyzed the flux at the origin a; = and found that the effect of the slow 
removal of the shutter is to suppress the DIT fringes in the current. It turns out that the 
Massey parameter (3 = jk/h identifies qualitatively two different regimes: for /? 1, one 
recovers the sudden approximation; but for /? ^ 5 the particle will experience a quasi- 
static barrier, no DIT-like transients in the current are observed, and the flux is indeed 
well-described by the equilibrium current 



AEt ~ h. 



(11) 



FWHM(£;) T > 2h. 



(12) 
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Fig. 4. Snapshots of the density profile of a chopped beam {h = m = 1) with 7 = 0.5 (a), 5 (b), and 20 
(c) taken in all cases Sit t = 2.50 with k = } . The solid line corresnonds to the sudden totally reflecting 
shutter (7 = 0), while tl 
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Fig. 5. Characteristic modulation of the density profile of a matter- wave beam of *^Rb atoms induced by 
a totally reflecting mirror which is displaced at constant velocity v > Vf^ {v^ = f .0 cm/s). The amplitude 
of the oscillations is maximized as v approaches the velocity of the beam V]^, when the density exceeds 
by more than 80% the classical case. In the limit ti ^ 00 the dynamics reduces to the diffraction in time 
setup with the lowest fringe visibility. 



J{x^O,t) = —T{E^—- 
m V 2m 



with 



T{E) 



E 



E^ 



(13) 



which is the transmission probability for a (5-barrier of strength ^/t. Nonetheless it is 
noteworthy that the dynamics of the density profile can be dramatically distorted with 
respect to the reference case (see Eq. Q) as shown in Fig.|4l due to the partial reflection 
from the shutter at short times. 

DIT is actually a delicate quantum phenomenon, easily suppressed or averaged out. 
An exception is the effect of an alternative removal of the shutter consisting of its dis- 
placement along the direction of the propagation of the beam. For a mirror moving with 
velocity v j42j . the constructive interference with the reflected components leads to an 
enhancement of DIT. In particular, the sharpness of the fringes in the density profile and 
their visibility are maximized whenever v — > pt/m, as illustrated in Fig. [5l 
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2.3. Dissipation 



It is clear that DIT requires quantum coherence. From this point of view it should 
come as no surprise that the presence of dissipation leads to its suppression. This was 
shown by Schuch and Moshinsky [38] using the phenomenological Caldirola-Kanai model 
[80|81j , which describes dissipation without explicitly considering the degrees of freedom 
of the environment, though it is roughly tantamount to the Caldeira-Leggett model in 
the Ohmic regime (82] ■ Within this approach, the fringes in the DIT analytical solution 
are washed out as the time of evolution goes by. 

2.4. The wave equations 

One may wonder if DIT could be observed in other systems, such as in the electro- 
magnetic field [5]. Let us consider the ordinary wave equation 



where the speed of light c. Moshinsky showed that no DIT phenomenon arises, the 
solution being 



with Lu = ck, which oscillates uniformly for t < x/c and vanishes elsewhere, in contrast to 
the DIT solution. Therefore, the DIT is absent in the free propagation of light. DIT may 
however be seen in constrained geometries, such as waveguides. Moreover, a relativistic 
approach to the diffraction in time using the Klein-Gordon equation shows that the 
probability density is restricted to the accessible region x < ct, at variance with the 
non-relativistic solution ip{x,t) = M{x,k,t) which is non-zero everywhere already for 
t = 0~^, see also 17.51 A similar discussion of DIT in the Dirac equation can be found in 
[83j . Finally, a mapping from space to time serves also to observe DIT from the vertical 
emission of a lasing mode in a laterally confined cavity [76j . More on this in Section 1101 

3. Transients in free expansion 

The free dynamics of a suddenly released matter-wave is frequently used to probe the 
properties of the initial state and its relaxation quantum dynamics [51 . The first studies 
generalizing the Moshinsky shutter setup in the context of neutron interferometry focused 
on rectangular wavepackets of the form e^''^X[o,L]: where the characteristic function in 
the interval [a, b] is X[a,b] = ©(^ ~ ~ ©((^ ~ EE]- We shall discuss the more realistic 
case of particle-in-a-box eigenstates. 



with n € N. Box-like traps have been implemented in the laboratory both in an all-optical 
way [49j and with microchip technology '50] . The expansion after suddenly switching off 



d^E{x,t) _ 1 d^E{x,t) 



(14) 



E{x,t) = e{t~x/c) e*(fc— *) 



1- 
2J' 



(15) 




(16) 
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Fig. 6. (Color online) Free expansion of the n = 1,2,4 eigenstates for a *^Rb atom released from a 
hard- wall trap of lO^im, see Eq. I I17II . 

the confining potential has been analyzed at short '9\ and arbitrary times [12] . Using the 
superposition principle and introducing fc„ = m:/L, the time evolved wavefunction is 



(j)„{x,t)^ / dx'Ko{x,t\x',t' ^0)^p{x',t' ^0) 



^^JtY. v[e"'''-''M(x~ L,vK,t) - M{x,vkr,,t)\. (17) 

Registration of the probability density at a given point as a function of time, allows 
to distinguish the following two regimes. Whenever the de Broglie wavelength is of the 
order of the box trap, which is to say n ~ 1, the particle expands following a bell-shaped 
profile with few transients as sidelobes. This is the Fraunhoffer limit of diffraction by a 
narrow slit in time. The DIT shows up for excited states n ^ 1, which corresponds to the 
so-called Fresnel regime, in which the effect of the two edges of the trap can be identified 
[9]. The space-time density profile is actually richer in structure as shown in Fig.[6l 

In what follows we shall compare the dynamics with the expansion from a harmonic 
trap. Let us first recall the spectrum of the harmonic oscillator (HO) trap, Sp{Hho) = 
{E^o ^ f^(^i ^ i/2)\l = 0, 1, 2, . . .}, and that of the hard-waU (HW) trap, Sp{Hhw) ^ 
{E^^ = {him/LY/{2m)\n = 1, 2, 3, . . .}. The condition for the n-th eigenstate of the 
HW trap to have the same energy than the {n — l)-th eigenstate of the HO leads to 

(2n-l)mL2- ^^^^ 
The eigenstates of the harmonic oscillator trap are well-known to be given by 

0„(x,t = O)= f — e- — i7„ f J— , (19) 

where Unix) are the Hermite polynomials [79 . Although the time-evolution of such 
eigenstates can be calculated for a broad variety of functional dependencies a;(t) [84], in 
the following we shall restrict ourselves to u>{t) — ujQ{~t) which is tantamount to the 
time-evolution under the free propagator Ko{x,t\x' ,0) in ([3]) for t > 0. The evolution is 
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Fig. 7. Density plots of tlie probability density for the dynamical evolution of a ^'^Rb atom in the tenth 
eigenstate released from a harmonic (left) and box (right) traps. L = 80 fim and w is chosen so that the 
energies of the eigenstate in both traps coincide. 

given in terms of a scaling law of coordinates up to an additional phase which will play 
no role in our discussion. Indeed, it reads |84| 

M^,t) = (pi = O) ("+V2)^ (20) 

where t :— dt' /b^{t') and b must be a solution of b+uj{t)b = with b{0) = 1 and 6(0) = 
0. An instantaneous turn-off of the harmonic potential at i = leads to b{t) = y/T+lvH^. 
Figure [7] shows the expansion of a single particle wavefunction suddenly released from a 
HO and a HW trap. If the particles are in the ground state of the trap the evolution is 
alike. However, for any excited state a remarkable difference appears. For the former case, 
the scaling law gives rise to a uniform spreading without affecting the initial quantum 
structure. As can be seen in the x-t density plot, the full with at half maximum (FWHM) 
becomes linear for any t ^ uj~^. By contrast, whenever the initial trap is box-like, the 
density profile exhibits a bifurcation into two main branches after the semiclassical time 

mL mL^ 

= ^ = (21) 

of expansion [12| . This is the result of the mapping of the underlying momentum dis- 
tribution to the density profile expected asymptotically. Indeed, the probability density 
in fc-space for a box eigenstate is |(^„(fc)p = ^p^[l — (— 1)" cos(fcL)]/(fc^ — fc^)^, and 

has two main peaks at A: = ±fc„ except for the ground state: |0i(A:)p, is a bell-shaped 
distribution centered at A; = 0. 

The more general dynamics in an expanding box has been described in 43 85j. The 
transient dynamics of particles released from a cylindrical trap has also been discussed 
in [55] and from spherical traps in [23110187] . We recall that the 3D s-wave case is tanta- 
mount to the ID problem in half space, with a hard-wall at a; = 0. One can exploit the 
method of images [16] to write down the propagator, 

X3^(r, t\r', 0) = Ko{r, t|r', 0) - Ko{~r, t\r' , 0), (22) 

and the wavefunction in the radial coordinate as 
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4>l''{r, t) - Mr, t) - M-r, t), r > 0, (23) 

where i^oir^t) = J^^dr'Ko{r,t\r' ,t' = 0)V'o(r',i' = 0), iJo{r',t' = 0) being the radial 
wavefunction of the initial state. The effect of interactions during the expansion will be 
considered in Section [5l 

3.1. Phase space representation 

3.1.1. The Wigner function 
The Wigner function [88] . 

oo 

W{x,p)-^^ I dy^{x + yrtP{x^y)e^'Py/\ (24) 

— oo 

is the best known quasi joint probability distribution for position and momentum. Nowa- 
days, it is possible to study it experimentally as has been demonstrated in a series of 
works |89I90I91I92I93I9^ . We recall that the marginals of the Wigner function are pre- 
cisely the momentum and coordinate probability densities, this is, J dxW{x,p) — \ip{p)\'^ 
and / dpW{x,p) = \'ip{x)\'^. 

In the context of the Moshinsky shutter, its time evolution was derived for the free 
case problem [5^ and in the presence of a linear potential |^ . The Wigner function for 
a cut-off plane wave initial condition isQ 

u/f 4- n^ 1 sm[-2x{pQ-p)/h] 

W(x,p;po,t = 0) ^ Q{-x), (25) 

TT PQ-p 

which clearly assumes negative values. If the potential is linear or quadratic, the time 
evolved Wigner function follows classical trajectories, 

W{x,p;t) = JJ dx.,dpiS[x ~ Xci{xt,pi,t)]6[p - pci{xi,pi,t)]W{xi,Pi;t = 0). (26) 

For (pS)) the result is limited to the classically accessible region of phase space. Moreover, 
we may take the limit ft ^ to find the classical distribution [29 96 9 71^5] . in which no 
hint about the quantum oscillations -hallmark feature of DIT- remains. 

Similarly, one can describe the time evolution of the eigenstates of a hard-wall trap. 
At the beginning of the experiment the Wigner function is given by |98] 

W4,„{x,p) ^Vn{x,p)xiO,L/2]{x) +Vn{L - X,p)x[L/2.L]{x), (27) 

where 

2 ( sr-^ sm[2(p/h+ i'mT/L)x] f 2mTx\ sin(2pa;//i) 1 

^" " ^ 1 4{p/h+ i^rnr/L) ''°%~^y 2p/ft J' 

Figure [5] shows the stretching of the Wigner function along the x-axis as time goes by. 
The momentum distribution, being bimodal for all n > 1, is responsible for the splitting 



Since 'ip(x,t = 0) = e^P°^^'^0(~x) is dimcnsionlcss, W(x,p) here has dimensions [p] 
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Fig. 8. Time evolution of the Wigner function for the fifth eigenstate of the box trap (L = lO^m, for a 
^'^Rb atom). As time goes by, the the quasi-probability distribution is stretched along the a;-axis. The 
gray scale varies from black to white as the function value increases, corresponding the background 
gray to zero value. Note in particular the presence of negative values (in black), telltale sign of the 
non-classicality of the state. 

in real space for t > f„ already noticeable at t — 5 ms. For a general Hamiltonian the 
propagation of a Wigner function obeys the quantum Liouville equation |88I99| 

dtW{x,p;t) + ^d,W{x,p;t) - ^ (2^ + 1)1 dl^+W{x)dl'+'W{x,p;t) = 0.(28) 

In such a way, the phase space description of transient effects can be extended to tunneling 
problems as well 1301100) . 

The Wigner function has been also proposed to study the quantum transient response 
to changes in an externally applied voltage in semiconductor structures [101] . see also 
Section [51 dynamics of polaron formation in compound semiconductors |102| . or for tran- 
sients in scattering processes [103] and classical-quantum comparisons [104] . 



3.1.2. Symplectic tomography 

A related representation in phase space is symplectic tomography [l05j . The wave 
function ip{x) or the density matrix p{x,x') can be mapped onto the standard positive 
distribution W{X, /z, u) of the random variable X depending on two real extra parame- 
ters, fi and v. The symplectic tomogram can be related to the Wigner function W{x,p) 



WiX,^l,ly)^ / W{x,p)S{X - nx- up) 



dxdp 
2'Kh 



(29) 



as its Radon transform. The parameters /i and v can be expressed in the form scos^, 
V — s~^sin0, so that s > is a real squeezing parameter and 6* is a rotation angle. 
Then the variable X is identical to the position measured in the new reference frame in 
the phase-space which has new scaled axes sz and s~^p and rotated by an angle 6. The 
Moshinsky solution in the tomographic representation reads [29] 



(30) 



where 
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^/2fJ.{flT + v) 

with T = ht/m. In the reference frame with {fJ-,v) = (1,0), the quantum tomogram 
reduces to the density profile in Eq. ([6]). 

3.2. Propagation of wavepackets with sharp boundaries. 

The sharp boundaries of the initial wavefunction considered so far, which also arise 
naturally in the Zeno dynamics under periodic spatial projections |106| , might cast some 
doubt on the realizability of DIT. The study of the generic time propagation of initially 
confined wave-packets has shed new light on the effects of the boundaries on quantum 
transients at short times [10711081109] . In the spirit of the Moshinsky shutter, let us 
consider a state iplxjt = 0) with support in the negative semiaxis {x < 0), and denote 
its Fourier transform by 

9ik) = J ^{x, t - Q)er^^^dx. (32) 

For any t > it follows that 



i^ix. t) = ^^-^ j g{k)M {x, k, t)dk. (33) 

An expansion of the Moshinsky function with respect to fc, allows to perform the integral 
in (|33p . and to express the solution in terms of derivatives of the initial state ip{x,t = 0) 
at the boundary x — 0, 



Expanding in a power series of {ht/2mY^^x — rj, one finds for ht/2mx^ <C 1, 



(34) 

v=o 



2TT'ip{x,t) ^ ^ -e"^ {t] - 2ijf{l + xdy) 

-Ari^W + xdy) + x^dl + . . ]}^{y, t = {})\y^^. (35) 

So, for short times, the value of the wavefunction depends exclusively on ^{x^ t — Q) and 
has a t^/^ dependence, whereas should ^{x,t = 0) vanish, higher order terms are to be 
considered. The effect of the exact wave packet boundary may however be washed out 
at longer times as demonstrated by the revival effect, see 14.41 

4. Beam chopping and formation of pulses 

4.1. Quantum transients in beam chopping 

Diffraction in time is generally associated with the release of a matter wave initially 
confined in a given region of space. Nonetheless, the opposite problem is of great phys- 
ical interest in intcrfcrometry, namely, the chopping of a beam by a shutter, or mirror. 
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Fig. 9. Snapshots of the density profile of a chopped beam {h = m = 1) with a) /c = at t = 5, b) fc = 1 
at t = 5, and c) A; = 1 at 4 = 20. The solid line corresponds to the totally reflecting shutter {u — > oo) in 
Eq. I I36I I. while the dashed line is for a finite <5-potential with i/ = 0.25, see Eq. I I38I I. 

Consider a monochromatic beam described by a travelling plane-wave ijj(x,t = 0) = 
gifcx g^^^ i-j-^g action of a totally reflecting mirror suddenly turned on at a; = and 
t = 0. For all t > the beam is effectively split into two different parts. Noticing that 
the corresponding propagator is that of a hard-wall at the origin, Kwauix,t\x' ,0) — 
Ko{x,t\x' ,0) — Ko{—x,t\x',0), the evolution of the beam is 



■0(X,<) = 



M{x,k,t) - M{-x,k,t), x<0 
M{-x, ~k, t) - M{x, -k, t), x>0, 



(36) 



where we have neglected the width of the mirror for simplicity. From the asymptotic 
properties of the Moshinsky function (see Appendix |B]) it is clear that the matter wave 

ipix) — 2ism{kx)e~'^^^ is eventually formed at the left of the mirror. The effect of 
finite reflectivity of the shutter can be taken into account considering a delta mirror 
V{x) = rj5{x) which represents a weak and narrow potential. The explicit dynamics is 
slightly more involved, for the propagator is |110| 



Ks{x,t\x' = Ko{x,t\x',Q) ^ xM{\x\ + \x'\,~iK,ht/ra), 



(37) 



with Kq{x, t\x' , 0) as in ([3]) and x = mij/h^. Using the integrals in Appendix iBl it turns 
out that 



ip{x,t) 



ikx—ihk t/2rti 



E 



ivk 



[M{\x\,vk,t) - M{\x\,-iK, t)], 



v=±\ 

which asymptotically, at long times, tends to 

fc 



'^(x,t) 



-ihk-'t/2m 



(38) 



(39) 



This is nothing but the stationary scattering state of the mirror with the energy of 



the initial beam (notice that 



and 



are the corresponding transmission and 
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reflection probability amplitudes, respectively). The dynamical buildup of this standing 
matter-wave after the action of the shutter is illustrated in Fig. [9l 

Further work on beam chopping and pulse formation has almost invariably been carried 
out using the "quantum source approach" , to which we now turn out attention. 

4.2. Formation of single pulses 

The quantum mechanical effects of beam chopping have been studied in time-dependent 
neutron optics |111I112| . The energy distribution of the initial beam is expected to be 
modified in this process in agreement with a time-energy uncertainty relation |23l24j . 
In this section we examine the formation of a single-pulse of duration r from a quasi- 
monochromatic quantum source iplx = 0, i) ~ g-^'^ot modulated according to a given 
aperture function x'-^''(i)[l] Here, the superscript denotes the creation of a single pulse 
|23I24I111I112I113I11I12IT5] . and the free dispersion relation holds, loq = hkl/2m. Many 
works have been carried out within this approximation in neutron interferometry [11111121114] , 
considering also a triangular |111I113| or more general apodized aperture functions [12], 
atom- wave diffraction [11], tunneling dynamics |58|115|116|117|33l34|32j . and absorbing 
media |36]. After the experimental realization of a guided atom laser |118j . the same 
formalism has been employed to describe the dynamics in such system [15]. 

For the sake of concreteness, let us focus on the family of aperture functions 

X«(t)-sin"(17t)e(t)e(r-t), (40) 

with f7 = 7r/r |12ll5j . The solution for an arbitrary aperture function is described in 
Section l43l We first note that the energy distribution varies for increasing n and T. The 
energy distribution of the associated wavefunction is proportional to (x = 0, w)p 

|119I15| . where 4>n\x = 0,w) = {2n)-'/^ J dt exp[i(a' — LOo)t)xn^ ■ For increasing n, the 
aperture function becomes smoother broadening the energy distribution. Similarly, the 
smaller the time for pulse formation r, the broader is its associated energy distribution. 
Indeed, whenever hujQT < 1, the distribution is shifted to higher energy components. 

As an example, let us consider the case Xo^''(^) = O(t)0(T — t) corresponding to a 
rectangular single-slit in time which is switched on and off at infinite velocity. The time 
evolution of the pulse is given by 

i{;l^^\x,ko,t;T) = [M{x,ko,t) + M{x,-ko,t)] 

-e(t - r)e^'''^"^[M(x, ko, t - t) + M{x, -ko, t - t)]. (41) 

li t < T, before the pulse has been fully formed, the problem reduces to that of a 
suddenly turned-on source, il)^\x,kQ,t < t) — M{x,kQ,t) + M{x,—ko,t), discussed in 
12.11 An explicit calculation of the energy distribution of the resulting pulse, taking the 
overlap of the wavefunction with the free particle eigenstates which vanish at the closed 
shutter [23124] . shows that 



^ Baute, Egusquiza and Muga 11191 have made explicit the relation between a source boundary condition 
where ■>p{x = 0, t) is specified for t > and the standard initial value problem in which t = 0) is 
specified. The simple "Kirchhoff boundary condition" ^(x = 0,t) = ©(t)e""o* is discussed further in l7.2l 
see also |120) . 
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V{E;Eq,t) =AfVE 



{E~E^f 

Af being a normalization constant. Such distribution is peaked at the energy of the initial 
beam Eq = Hujq if Eqt > h and shifted to higher energies otherwise H^. Indeed, for 
some measure of the width AE it leads to the time-energy uncertainty relation AEt ~ h 
[23l24l25ll2j . 

It is illuminating to consider a slightly more complicated aperture function, such as 
the sine-square (Hanning) function, 



1 - - cos(2ffi) 



e{t)e{T-t). 



(42) 



Defining fc^ — ^ 2m{uJo ± 2f2) / h with /3 = ±1, the resulting pulse is described by the 
time-dependent wavefunction 



tp2\x,ko,t;T) 



/3=±1 



(43) 



where the effect of the apodization is to subtract to the pulse with the source momentum 
two other matter- wave trains associated with kp, all of them formed with rectangular 
aperture functions. The effect of creating the pulse by a smooth modulation as opposed 
to the sudden rectangular aperture function, is to suppress the diffraction sidelobes in 
the density profile, at the expense of broadening the corresponding energy distribution. 
These two features are the key elements of apodization in classical optics and Fourier 
analysis, and hence the effect was dubbed apodization in time [12ll5j . 

4.3. Arbitrary modulation of a matter-wave source 



Though the dynamics of matter-wave sources associated with different aperture func- 
tions has been worked out in detail [2411 1111 1211 1311 lll2j . it is desirable to tackle the 
general case. In order to do so, we consider an arbitrary aperture function x{t) which is 
zero outside the interval [0, r] [15j . An analytical result can be obtained in such case by 
means of Fourier series. Noting that 



xit)= E 



-e{t)Q{T-t) with 



Cr = 



'dt, 



(44) 



the wavefunction at the origin can be written as the coherent superposition of its Fourier 
components modulated with a rectangular aperture function 



iP^(x :^0,ko,t)^ E Cre~"^'-*Q{t)e{T -t), 

r—~oo 

with = ^0 + 27rr/r. Clearly, the subsequent dynamics for x, i > reads 

oo 

ip^{x,ko,t;T)= E Cripo\x,kr,t;T), 



(45) 



(46) 



r— — OO 
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Fig. 10. Apodization in time. Removal of the sudden approximation in the switch of a matter wave source 
leads to the suppression of oscillations in the beam profile characteristic of the diffraction in time. The 
^'^Rb beam profile is shown at t = 20ms, and hko/m = Icm/s. For increasing t the amplitude of the 
oscillations diminishes and the signal is delayed. 

where ip'ip (x,k,t]T) is defined in Eq. (|4ip . and hk^ = -v/^toHov, with the branch cut 
taken along the negative imaginary axis of ujr- The upshot is that the dynamics of a 
source with an arbitrary aperture function can be described in terms of the coherent 
superposition of rectangular pulses associated with each of the Fourier components. 

4.4. Diffraction-in-time revival after slow switching 

The dynamics of a matter-wave source can be tamed by switching it slowly. For the 
sake of concreteness, we introduce the continuous switching function 

(-0 if t < , 

^sn,^tch^^■^^] ^(^)^ if 0<t<r, (47) 

U, if i>T, 

with a characteristic time scale r, which we shall refer to as the switching time. The 
family for which — sin" ilgt with n = 0,1,2 and f2s — 7r/2T was considered in 
[T5.121j . where it was shown that for fixed r, the apodization of the beam increases with 
the smoothness of this is, with n. Note that for n = 0, one recovers the sudden 
aperture Xo(0 — Q(^) which maximizes the diffraction-in-time fringes. We shall consider 
a sine-square modulation in what follows (n = 2), and concentrate on the dependence on 

T. 

Increasing the switching time r leads to an apodization of the oscillatory pattern. As 
a result the fringes of the beam profile are washed out, see Fig. [TOl In this sense, the 
effect is tantamount to that of a finite band- width source [34j , and as already remarked, 
this is the essence of apodization (in time), the suppression of diffraction by means of 
broadening the energy distribution [122ll9j . 

However, the apodization lasts only for a finite time after which a revival of the diffrac- 
tion in time occurs, as shown in Fig. 1111 The intuitive explanation is that the intensity of 
the signal from the apodization "cap" associated with the switching process during the 
interval [0, r], decays with time, whereas the intensity of the main signal (coming from 
the step excitation for t > t) remains constant. For sufhciently large times, the main sig- 
nal, carrying its diffraction-in-time phenomenon, overwhelms the effect of the small cap. 
Indeed, from the linearity of the Schrodinger equation, it follows that the wavefunction 
is the sum of a "half-pulse" term associated with the switch released during the interval 
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Fig. 11. Revival of the diffraction in time. During the time evolution of an apodized source, there is 
a revival of the diffraction in time in the time scale tr = ujot'^ {t = 1ms, Hko/m = 0.5cm/s, and 
tr = 16.8ms for a *^Rb source which is switched on following a sine-squarc function). 

[0, t] and a semi-infinite beam suddenly turned on at time t = t. The revival time can 
be estimated to be |15|121) 

U^ojot'^, (48) 

and plays a similar role to the Rayleigh distance of classical diffraction theory, but in 
the time domain jl23j . The smoothing effect of the apodizing cap cannot hold for times 
much longer than tj.. This is shown in Fig. where the initially apodized beam profile 
eventually develops spatial fringes, which reach the maximum value associated with the 
sudden switching at i « 2tj.. 

4.5. Multiple pulses 

4.5.1. Application to atom lasers 

Knowledge of the dynamics associated with a single-pulse if^''^'' {x,t) can be used to 
tackle matter wave sources periodically modulated in time. Frank and Nosov first de- 
scribed the periodic action of a quantum chopper in ultracold neutron interferometry 
leading to multiple rectangular aperture functions [13|14j . More recently, del Campo, 
Muga and Moshinsky [15] found the general dynamics of an atom source under an ar- 
bitrary aperture function as described in 14.31 ^^id applied it to the description of an 
atom laser in the non-interaction limit. The experimental prescription for an atom laser 
depends on the "out-coupling" mechanism. Here we shall focus on the scheme described 
in |124ll25j . The pulses are extracted from a condensate trapped in a magneto-optical 
trap (MOT), and a well-defined momentum is imparted on each of them at the instant of 
their creation through a stimulated Raman process. The result is that each of the pulses 
'0^^'' does not have a memory phase, the wavefunction describing the coherent atom laser 
being then 

N-l 

<l)^^\x,ko,t) = J2 i'^'H^,ko,t~jT;r)e{t-jT). (49) 
j=o 

Such modulation describes the formation of N consecutive x^^^-pulses, each of them 
of duration r and with an "emission rate" (number of pulses per unit time) l/T. The 
Heaviside function 8(t — jT) implies that the j-th pulse starts to emerge only after jT. 
We can impose a relation between the out-coupling period T and the kinetic energy 
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Fig. 12. Dynamics of the density profile of an atom laser a) coherently constructive, b) incoherent case. 
Note that Xd = hkot/m is the classical trajectory. In all cases the norm is relative to the incoherent 
case, as explained in the text. (^"^Na source, modulated by a sine-square function with r = 0.833ms, 
T = 0.05ms, moving at 3.73cm/s, for which the revival time is 174.4 ms.) 

imparted to each pulse huja, such that woT — In: if I is chosen an even (odd) integer the 
interference in the beam profile is constructive (destructive). 

The diffraction in time is a coherent quantum dynamical effect, so we may expect 
it to be affected by noise in the environment, ff the phases between different pulses 
are allowed to fluctuate, the interference pattern is blurred [126] . Moreover, if there is 
no phase coherence between different pulses, the resulting density profile becomes the 
incoherent sum of the single-pulse densities, p{x, t) = \ip''^^x, ko,t—jT; T)\'^Q{t— jT), 
which we shall use as a reference case. 

One advantage of employing stimulated Raman pulses is that any desired fraction of 
atoms can be extracted from the condensate. Let us define s as the number of atoms 
out-coupled in a single pulse, The total norm for the A^-pulse incoherent atom laser 
is simply Mine ~ sN. A remarkable fact for coherent sources is that the total number 
of atoms A/'c out-coupled from the Bose-Einstein condensate reservoir depends on the 
nature of the interference [127] . Therefore, in Fig. [T^] we shall consider the signal relative 
to the incoherent case, namely, {x,t)\'^ /Afinc- 

Note the difference in the time evolution between coherently constructive and incoher- 
ent beams in Fig. [121 At short times both present a desirable saturation in the probability 
density of the beam. However, in the former case a revival of the diffraction in time phe- 
nomena takes place in a time scale loot'^, similar to the smooth switching considered in 
14.41 whereas the latter develops a bell-shape profile. This revival time can actually be re- 
lated to the frequency of the BEG reservoir trap uJtrap, which determines the width of the 
out-coupled pulses [121] . Assume that the atoms are out-coupled from the n-th longitudi- 
nal mode of the reservoir trap, whose spatial width is ^„ = [{n + l/2)h/mujtrap]^^^ ■ Using 
ujo — ?ifc§/2m, and the semiclassical relation t„ = m^„/?ifco, the revival time becomes 




(50) 
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For a non-interacting BEC, one can take n — 0. Therefore, when the trap of the BEC 
reservoir is very tight, DIT will arise ever since the creation of the atom laser beam, while 
if the trap is wide, so will be the out-coupled pulses and apodization will suppress the 
oscillations on the density profile for times smaller than the revival time. Conversely, if 
the initial front of the beam profile is known to be smoothed out in a given length scale 
^, it can be shown that density modulations arise in the time scale m^^/fi,^. 



4.5.2. Coherent control 

As an outlook, the techniques of "coherent control" |128|129ll30j . could merge with the 
atom laser to design matter-wave pulses for specific aims. A remarkable example is the 
deterministic generation of arbitrary de Broglie wave- fronts in atom lithography [131j . 
Coherent control has been applied in recent years to a wide range of systems in atomic, 
molecular, solid-state physics, semiconductor devices, biology or quantum information. 
The basic idea is to make use of quantum interference to manipulate the dynamics with 
the aid of a control field, typically formed by laser light. Examples of objectives so far are 
selective bond breaking, choosing a reaction pathway, or performing a quantum infor- 
mation operation. The experimental output may be used in the optimization procedure 
with closed- loop learning control techniques. Instead of a control field made of light, we 
envision the formation of control matter-fields by combining matter-wave pulses with 
tunable intensity, energy, timing and phase. An elementary example is the formation of 
a continuous atom laser by constructive interference, as discussed above, and other ap- 
plications may follow in which the elementary pulses and their transients will be helpful 
building blocks for analyzing the resulting matter field. 



4.5.3. Ultra-short laser pulses 

Ultrashort laser pulses in femtosecond or attosecond time domains are a basic tool in 
coherent control research and applications. These nowadays manipulable pulses are able 
to induce many transient phenomena for nuclear or electronic motion. Using pump and 
probe or streaking methods, the transient dynamics following a system excitation may 
be recorded and/or steered. A detailed account is out of the scope of this review and 
impossible to summarize here, but we shall mention briefly some recent examples to show 
that there exists a link with work on elementary quantum transients which has not been 
fully exploited and is worth pursuing. 

Strong field ionization of atoms via tunnelling can be nowadays followed with attosec- 
ond resolution so that some aspects of electron tunnelling (see others in Section [7|) can 
be observed |132ll33j . Impressive "streaking techniques" are used for this purpose, in 
which emission times and momenta of electrons are mapped, thanks to the action of the 
short pulses on the ejected electron, into final ion and electronic momenta which are 
readily observable with time-of-flight detectors. The emission may also be controlled so 
as to produce double or multiple time slits with few-cycle pulses and the correspond- 
ing interference-in-time phenomena |134| . Simple models for apodized pulses or time- 
dependent shutters (see e.g. Section 12. 2p could be contrasted with these experimental 
procedures and may help to interpret the observed results and design new experiments. 
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5. Dynamics of ultracold gases in tight-waveguides 



The models described so far do not explicitly consider quantum statistics and rig- 
orously only apply to single-particles. Therefore, even though the use of semi-infinite 
beams and quantum sources is suitable to describe neutron optics and ideal atom waves, 
a more rigorous approach is desirable. As already pointed out by Stevens, a consistent 
interpretation for several -non-interacting- particles requires the appropriate quantum 
symmetrization |135j . Moreover, there is nowadays a flurry of theoretical and experimen- 
tal work dealing with ultracold atomic gases in the non-linear and strongly interacting 
regimes. 

Ultracold atoms in waveguides, tight enough so that the transverse degrees of freedom 
are frozen out |136|137j . are accurately described by the Lieb-Liniger (LL) model, 

"~2^E5^+^'i^ E ^i^^-^i)^ (51) 

1=1 * l<i<j<N 

exactly solvable by coordinate Bethe ansatz [138j . This model describes N bosons in 
one-dimension interacting with each other through a short-range (5-function. The use of 
such pseudo-potential is justified provided that the scattering is restricted to s-wave type 
in the zero-energy limit of relevance to ultracold temperatures, when the details of the 
true interparticle potential become unimportant [139]. 

The cloud can then be characterized by the interaction parameter 7 — 'mgijjL/h?N, 
where giD is the one-dimensional coupling strength, L the size of the system, m and N 
the mass and number of atoms respectively; 7 can be varied 1136 allowing to realize the 
mean- field regime (7 <C 1) or the Tonks-Girardeau regime (7^1) _140j . 



5.1. Dynamics in the Tonks-Girardeau regime 



We shall next consider the dynamics of strongly interacting bosons confined in a tight- 
waveguide, the so-called Tonks-Girardeau (TG) gas. In the TG regime, bosons mimic 
an effective Pauli exclusion principle as a result of the repulsive interactions rather than 
the quantum statistics itself. The Fermi-Bose (FB) duality |141I142I143I144| provides 
the many-body wavefunction of N strongly interacting bosons from that of a free Fermi 
gas with all spins frozen in the same direction. The Fermi wavefunction, antisymmetric 
under permutation of particles, is built as a Slater determinant, with one particle in each 
eigenstate (/>„ of the trap, 

i;F{xi,...,XN) = --^det'^f,^j^(t>n{xk). (52) 

Whenever the position of two particles coincide, the wavefunction ipp vanishes as a 
consequence of the Pauli exclusion principle. This contact condition is the desirable one 
in the TG gas due to the effectively infinite repulsive interactions. Symmetrization can 
then be carried out "by hand" , applying the so-called "antisymmetric unit function" 

A= Y]. s,gn{xk-Xj), (53) 

i<j<k<N 
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so that 



iPb{xi, ■ ■ . ,XAr) = A{xi, . . . ,xn)^f{xi, ■ ■ ■ ,xn)- (54) 

This is the Bose-Fermi map introduced by Girardeau in 1960 |141| . A remarkable advan- 
tage of this mapping is that it holds for time dependent processes (governed by one-body 
external potentials), since the A operator does not include time explicitly [144|142j . The 
glaring upshot is that as far as local correlation functions are concerned, to deal with the 
many-body TG gas it suffices to work out the single particle problem, since 

\iPb{xi, . . . ,XN;t)\'^ = \ipFixi, ■ . .,XN;t)\^- (55) 

In particular, from the involutivity of the A operator {A'^ = 1) and the fact that 
{(f>n\U''U\<j)„i) = Snnn where U is the time-evolution operator, it follows that the time- 
dependent density profile can be calculated as |144j 

N 

p{x,t) ^ N / \tpB{x,X2,.-.,XN;t)\'^dX2 - ■ -dXN ^^^\4'n{x,t)\'^. (56) 

Once the single-particle time evolution is known, we are ready to study the spreading of 
the TG gas through its density profile, Eq. (|56p . 

For the HO trap, several studies have been carried out describing the scaling law 
governing the expansion of the TG gas |54I145| as well as the dynamical fermionization of 
the system exhibited in the momentum distribution [1461147] . Here, we shall focus on the 
transition to the ballistic regime in the expansion from both HO and HW traps. Moreover, 
we shall be interested in low number of particles N, where the mean-field approach is not 
accurate, missing the spatial anti-bunching in the density profile |148|144j . If we are to 
compare the expansion from both traps, it seems natural to consider that, for the same 
number of particles, the total energy of the TG gas is to be the same. Then, the following 
relation must hold between the harmonic frequency and the length of the box, 

_ h7r^N+l){2N + l) 

- 6^^^ ■ ^^^^ 

For the HO trap, the similarity transformation 

p{x,t) = ^=J==p (^=^==,0) (58) 



has been shown to hold, where the ballistic regime sets for t ^ uj~^ [51]. However, 
as we have already discussed, no such simple expression can be obtained for the hard- 
wall trap, exhibiting the dynamics a lack of self-similarity (which should be clear once 
the transient features entailed in Eq. p7|) are recognized), see Fig. [131 The ±pjv = 
±hNT:/ L components govern in this case the width of the expanding cloud for t >tN — 
mL^/2NTrh. 

Figure [13] shows the variation in time of the full width at half maximum (FWHM) for 
a TG gas expanding from both types of traps. The transition to the ballistic regime is 
sharp for the HW trap a,t t^, whereas for the HO happens gradually, and only sets for 
t ^ uj^^. This fact, together with transient versus self-similar dynamics, points out the 
relevance of the confining geometry. 
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Fig. 13. Density plots of the probability density p for a TG gas composed of Af = 10 atoms of °'Rb 
released from a harmonic (left) and box-like (right) trap. L = 80 ^m, and the frequency w is chosen 
according to Eq. I I57I I. 
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Fig. 14. Time dependence of the FWHM of a cloud of A'^ = 10 *^Rb atoms in the Tonks-Girardeau 
regime, released from a box-like (L = 80 fim, solid line) and harmonic (dashed line) traps. Note that the 
transition to the ballistic regime is sudden when the gas is released from a box-like trap, but smooth if 
initially confined in a harmonic trap (parameters as in Fig. I13I I. 

The density profile we have described so far in the TG regime, is actuahy the same that 
for spin-polarized non-interacting fermions. We note that for the latter the momentum 
distribution 

n(fc) = ^ / dxdye'''^'=-y'>p{x,y) (59) 

remains invariant under time evolution provided that the momentum operator commutes 
with the purely kinetic Hamiltonian. Since the reduced single-particle density matrix 
(RSPDM) of spin-polarized fermions is 

p{x,y) = N / dX2 ■ ■ ■dXN'tpF{x,X2, ■ ■ ■ ,XN)*1pF{y,X2, ■ ■ ■ ,xn), 
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N 



= ^C(x)0„(y), (60) 

n=l 

the momentum distribution is the sum 

N 

nF(fc) = ^ |^„(A:)|2, (61) 

(j>n denoting the Fourier transform of However, the momentum distribution of 

an expanding Tonks-Girardeau gas exhibits a transient dynamics. An efficient way to 
compute the RSPDM has been recently derived by Pezer and Buljan |149j expHcitly 
using the Bose- Fermi map and the Laplace expansion of the determinant, see also [150) . 
Thanks to the orthonormality of the single-particle eigenstates, the RSPDM is given by 

N 

pTG{x,y)^ ^ (l)i{x)Ain{x,y)(j>n{y), (62) 

l,n=l 

where 

A(a;,y) = (p-i)^detP, (63) 

and the elements of the matrix P are 

y 

Pin ^ 6i„ - 2 f dz(j)Uz)Mz) (64) 



for X < y with no loss of generality. A generalization of this result for hard-core anyons 
can be found in |151| . 

From ((59|) and (|62|) . the time-evolving n{k,t) can then be found. For t = it presents 
a prominent peak a,t k — 0, whereas as the cloud expands and the particles cease to 
interact, the interaction energy is gradually transformed into kinetic energy, and the 
asymptotic momentum distribution is the quasi-momentum distribution, namely, the 
flat one typical of the dual Fermi system. This is the so-called dynamical fermionization 
[146|147j . illustrated in Fig. [15] for a TG gas expanding from a box-like trap. 

In conclusion, strongly repulsive interactions tend to suppress quantum transients in 
the density profile, while enriching the dynamical picture in the momentum space. Con- 
versely, attractive interactions have been shown to lead to an enhancement of density 
modulations in both expanding clouds and propagating beams |152j . 



5.2. Finite interactions: dynamics of Lieb-Liniger gases 



Whenever the interactions are finite, the exact dynamics of ultracold gases becomes a 
highly non-trivial subject. The first steps to cope with it made use of a hydrodynamical 
approach |137I54I145| and numerical simulations, but only recently have the required 

^ Note the errata of the pubhshcd version (Physics Reports 476, 1-50, 2009), where in Eq. (60) "0* 
should read and in Eqs. (62) and (64) i/>* should read </)*, as corrected here. 
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Fig. 15. Dynamical fermionization of TV = 5 ° Rb atoms in the Tonks-Girardcau regime, released from 
a box (L = 80 fim). The asymptotic momentum distribution is that of a spin-polarized ideal Fermi gas. 

analytical techniques been developed |153I154I155| . Within the Lieb-Liniger model the 
wavefunction of the gas ■0_b(2^1j • ■ • describes N identical (5-interacting bosons in 

one-dimension governed by the Hamiltonian ([5T|) |138j . It is customary to express the 
corresponding time-dependent and many-body Schrodinger equation (TDMBSE) with 
"Lieb-Linigcr" units (2m = 1, ?i = 1) in the form 



We assume that the gas is initially confined in an external potential V{x)^ which is 
suddenly turned off at t = 0, so that the cloud starts expanding. Therefore for i > the 
domain of each coordinate is the whole real line Xj S M. 

The wavefunction V'b is totally symmetric under permutations of particles, and hence 
it suffices to focus on the fundamental sector Ri : xi < X2 < ■■■ < xn, where the 
TDMBSE reduces to 



The Lieb-Linger contact conditions are responsible for creating a cusp in the wave func- 
tion when two particles touch. This can be expressed as a boundary condition at the 
borders of i?i [T38] . 



These boundary conditions can easily be rewritten for any permutation sector. In the 
TG limit (c — > oo), the cusp condition implies an effective exclusion principle whereby 
the wavefunction identically vanishes when two particles are in contact [141j . Exact 
solutions of the time-dependent Schrodinger equation (I65|) can be obtained by using 
a Fermi-Bose mapping operator [15611541155] acting on fermionic wave functions. As- 
sume 4'f{xi, . . . ,XN,t) is an antisymmetric (fermionic) wave function, which obeys the 




(65) 




(66) 




(67) 
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Schrodinger equation for a non- interacting Fermi gas, i-^f- = ^ X]i=i qx'-^ • Then, the 
wavefunction of the Lieb-Liniger gas can be written as 

^PB^Kda4>F, (68) 

using the generahzed Fermi-Bose mapping operator 

1 



n 

i<»<i<JV 



(- 


--)] 




dxi J 



(69) 



up to the normahzation constant A/'c [156J. This surprising result states that the time evo- 
lution of a Lieb-Liniger gas at a given time t can be computed from that of non- interacting 
fermions at the same time, and only then take into account the interactions through the 
generalized Fermi-Bose map. The applicability of this map is not limited to the expan- 
sion dynamics of an excited Bose gas initially confined in an external potential V{x), for 
it actually describes the ground state whenever [Oc, V{x)\ « 0. Nonetheless, as the gas 
expands in ID, the linear density decreases while the effective coupling constant becomes 
larger, prompting the system to enter the strongly interacting regime. Indeed, the asymp- 
totic form of the wavefunction can be probed to possess the Tonks-Girardeau structure 
[155j . However, the properties of this asymptotic state can considerably differ from the 
properties of the Tonks-Girardeau gas in the ground state of an external potential [155j . 
We note that such dynamics is restricted to free expansion, while the inclusion of specific 
boundary conditions (periodic, Dirichlet, Neumann) and its extension to other quantum 
statistics (interacting fermionic or anyonic systems) and multi-component, multi-channel 
problems remain as interesting open problems. 

5.3. Finite interactions: mean-field approach 

By now, it should be clear that suddenly released matter-waves initially confined in a 
compact support exhibit quantum transients. Nonetheless, the effect of strongly repulsive 
interactions tends to smooth out the density profile in the resulting dynamics, as in the 
Tonks-Girardeau gas. In this section we discuss the effect of the interactions on quantum 
transients within the mean-field approach. 

Weakly interacting ultracold gases in ID are then described by the Gross-Pitaevskii 
equation, as a result of approximating the wavefunction of the Bosc-Einstcin condensate 
by the Hartree-Fock ansatz '^{xi, . . . ,X]s) = Y\f=i4>(.^i)- For an effectively ID Bose- 
Einstein condensate, in the presence of some external potential V{x) (i.e. the trap) the 
condensate wavefunction <I>(a;) = \/N<j){x) obeys the ID time-dependent Gross-Pitaevskii 
equation 



dt 2m dx"^ 



V{x)+giDmxt 



$(a:), (70) 



where m is the atomic mass, and giD — ~2h^ /raauj the effective ID coupling parameter, 
aiD being a known function of the three-dimensional scattering length 1136'. In the mean- 
field regime, for low enough temperatures, the phase fluctuations can be suppressed [157j . 
Making use of time-dependent shutters a way to enhance the diffraction in time was shown 
in [42j . An alternative way is to make the interactions in the system attractive, (say, by 
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means of a Feschbach resonance [51j ) as discussed in |152| . Quantum transients arising in 
suddenly released matter-waves can then be enhanced in systems governed by attractive 
interactions and suppressed in the repulsive case. 

Indeed, it is worthy to consider the case in which the mean-field interaction domi- 
nates over the kinetic energy (but at boundaries). In such a case, the so-called Thomas- 
Fermi approximation assumes that the kinetic energy can be neglected in the Hamilto- 
nian so that the time-independent Gross-Pitaevskii equation reads ii^tfIx) — [V{x) + 
giD\^TF{x)\^]^TF{x), where ^ is the chemical potential. The Thomas- Fermi wavefunc- 
tion is then given by ^tf{x) = [(/i — V{x))/ gio]^^'^ whenever ^ > V{x) and zero else- 
where. The dynamics of the Thomas-Fermi density profile nxFix) = |$TF(a;,t)P after 
switching off the trap obeys well-known scaling laws with no transient behavior whatso- 
ever |158I159| . For instance, for the case of a harmonic trap V{x,t) = mw(t)^x^/2. 



(71) 



where the scaling coefficient satisfies the differential equation b = uP' /iP' subjected to the 
initial conditions 6(0) = 1, and 5(0) = 0. Clearly, the dynamics entailed in Eq. iffTjl lacks 
any DIT-related transient. 



5.4. Turning interactions on and off in a Bose-Einstein Condensate 

The world of cold atoms provides many interesting transients due to the possibility to 
turn on and off interatomic and/or external interactions. Ruschhaupt et al. [160] have 
examined the short-time behavior of a Bose-Einstein condensate when the interatomic 
interaction is negligible for the preparation in the harmonic trap, and strongly increased 
when the potential trapping is removed. A quantum interference effect in momentum 
space is then found in the Thomas-Fermi regime: the momentum distribution expands 
due to the release of mean field energy and the number of peaks increases with time one 
by one because of the interference of two positions in coordinate space contributing to the 
same momentum. The effect is stable in a parameter range and could be observed with 
current technology. Interestingly enough, similar abrupt changes in time of the coupling 
constant giD [161] and density perturbations [162j can also induce self-modulations in 
the density profile and shock waves. The appearance of shock waves is a phenomenon 
also present in classical non-linear equations and is out of the scope of this review. We 
refer the reader to the good existing works on the topic [16311641165] . More generally, the 
dynamics after a quantum quench in the Hamiltonian in a many-body system exhibits 
non-trivial quantum transients, i.e. see [166] and reference therein. Nonetheless, one has 
to bear in mind that the nature of this type of transient is dramatically different of those 
related to the DIT, for the non-linear interactions entail a dispersion relation different 
from the free one. 



5.5. State reconstruction 



The advances in the dynamical description of ultracold atoms in tight-waveguides have 
paved the way to the study of transients of Lieb-Liniger gases [153|154|155] . Nevertheless, 
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the sudden quench of interactions using Feschbach resonances |51) combined with a switch 
of the confining potential stands for its apphcations in tomography of trapped ultracold 
gases [55j . The evolution of the essentially free single-particle density profile, no(x,t), 
allows us to obtain the reduced density matrix at t = 0, 

p{k, k')^— [ fio(k' -kMk' ~ kle'^^^^^^^dt, (72) 
m J 

where the Fourier transform of the density profile reads 

Mk, i)^^ J Mx, t)e-'^''dx. (73) 

Alternatively, the RSPDM can be diagonahzed [167j as p{x, x' ]t = 0) = Y^- Xjip* {x)(pj {x') 
in terms of the orthonormal natural orbitals ipj (x) with occupation numbers Xj > sat- 
isfying J p{x,x')(pj{x)dx — Xjipj{x') and — Under free evolution, having set 
up c{t > 0) = 0, the density profile reads 

no(x,t) = p{x,x;t) ^ ^Xj\ipj{x,t)\^ , (74) 



with ifj{x,t) = (27r)-i/2 J ^f^~.(^i^-^^tkx-ihk''t/2m^ which, using Eq. ([72]), leads to the den- 
sity matrix p{k,k') — Xjip*{k)tpj{k'). Experimental measurements are restricted to 
t > 0, limiting the possibility of state reconstruction by means of Eq. ([7^ . but assuming 
time-reversal invariance or using spatial symmetries of the system, the integral in time 
can be extended over the whole real line [22] ■ More generally, once the interactions in 
the system are negligible, the reconstruction of the RSPDM of the initial quantum state 
from the cloud dynamics in a potential V{x) has been successfully addressed by Leon- 
hardt and coworkers within the optical tomography |168ll69j . Remarkably, the Wigner 
function of non-interacting Helium atoms has been experimentally measured by looking 
at the time-evolution of the density profile [92], from which higher order correlations of 
the trapped gas can also be inferred |170j . 



6. Transient effects with external potentials 

Transient effects have been considered so far mostly for motion in free space but new 
and interesting features arise in the presence of external potentials. Most work has been 
carried out within the sudden approximation for the shutter, namely, its instantaneous 
release. The quantum transients with short-range external potentials often admit an 
exact solution, as is the case for a (5-barrier |171I110>172. 173.174.30. 67 108J. the time- 
dependent (5-barrier [6140] . and the Kroning-Penney lattice used to mimic the behavior 
of electrons in crystals [175^ Similarly, it is possible to tackle multichannel problems 
such as the excitation dynamics of atomic wavepackets interacting with narrow laser 
beams |17) . Much work on these potential-dependent transients has been carried out to 
understand time dependent aspects of tunneling and will be reviewed in Section [T] 
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6.1. Extended quasi-monochromatic initial states 



Here we review an approach, first considered by Garcfa-Calderon and Rubio I61j, that 
leads to an exact analytical solution to the time-dependent Schrodinger equation with 
cutoff wave initial conditions for arbitrary potentials of finite range. The approach in- 
volves the complex poles and residues (resonant states) of the outgoing Green's function 
of the problem and renders the time-dependent solution both along the internal and the 
transmitted regions of the potential. These authors solve the time-dependent Schrodinger 
equation 

ih-^^-H^^/j{x,t)^0, (75) 

where H = —{h^/2m)d^/dx'^ + V{x), and V{x) describes a potential of arbitrary shape 
extending from x ~ to x — L, with the plane wave initial condition 

ip{x,t = 0) = e'''^'ei-x), (76) 

which is the same initial condition considered by Moshinsky for the free case, see Eq. 
([H). Garcia-Calderon and Rubio [ST] used Laplace transform techniques to obtain along 
the internal region of the potential the expression [61| 

ip{x, k, t) = (i){x, k)M (0, fc, t) - 

OO 

^ </.„(a;)M(0,fc„,i); (0 < x < L), (77) 

n— — oo 

where (j){x,k) refers to the stationary solution, and 4>n{x) = iun{0)un{x) / {k — fc„) is 
given in terms of the poles {kn} and resonant states {u„(a;)} of the problem, discussed 
in Appendix [E) Notice that along the internal region of the potential, the Moshinsky 
function does not depend on x and hence it does not correspond to a propagating solution. 
Similarly, along the transmitted region, the above authors provided the solution 

oo 

%l}{x,k,t) = T{k)M{x,k,t)- ^ TnM{x,kn,t); (x > L), (78) 

n— — oo 

where T{k) is the transmission amplitude and T„ = m„(0)M„(L)e^*'^"^/(A: — fc„). Fur- 
ther notice that in the above two equations the function M{x, fc„, t) corresponds to the 
Moshinsky function for the complex value k — kn- 

In the absence of a potential, the solution given by Eq. (|77|) vanishes exactly and in 
Eq. ([75)l . T{k) — 1 and the resonant sum vanishes as well, so that the solution becomes 
the free solution. In [^1], it is shown that the exact solutions given by Eqs. ([77]) and ([75]) 
satisfy the corresponding initial conditions, i.e., they vanish exactly for t = Q. Similarly, 
it is also shown that at asymptotically long times the terms M{x,kn,t) that appear in 
the above equations, tend to a vanishing value, while M{x,k,t) tends to the stationary 
solution, as first shown by Moshinsky [5]. Along the internal region 

ip{x, t) - (/)(x, k)ejip{-iEt/h), (79) 
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Fig. 16. Plot of the building up of \^{x,t)\'^ along the internal region of a double-barrier potential for 



several times with the parameters given in the text. At t 
from the stationary solution. 

and along the external (transmitted) region, 
ip{x,t) ^ T(k)exp{ikx)exp{—iEt/h). 



15 ps, \ip(x,t)\ becomes indistinguishable 



(80) 



The above solutions allow us to study the dynamics along the full time span from zero to 
infinity. In [6T] these solutions were used to examine the transient behavior of a double 
barrier resonant structure possessing an isolated resonance level and an initial plane wave 
with energy equal to the resonance energy of the level. The calculations showed that for 
times roughly of a fraction of a lifetime onwards the single-resonance approximation to 
the time-dependent solutions given above was in good quantitative agreement. Figure [TCI 
exhibits the building up of the probability density ["(/"(a^, i)P along the internal region of 
the double-barrier system (DB) for several times until it reaches the stationary solution. 
The parameters of the DB system, are typical of semiconductor resonant structures [207j : 
barrier heights V — 0.23 eV, barrier widths b — 5.0 nm, well width w — 5.0 nm and 
effective electron mass m — 0.067me, with rUe the electron mass. The isolated resonance 
level is characterized by the resonance parameters £i — 0.08 eV and Fi = 1.0278 meV. 
Figure [17] provides a comparison of the corresponding transmitted probability density 
with the free evolving solution at long times, which exhibits a delay time of the order of 
the theoretical estimate r = 2?i/Fi. 

The fact that the plane wave initial state jumps from unity to zero at a; = implies 
that the resonant sums in Eqs. ([77]) and ([75]) have a slow convergence. For this reason, 
most subsequent work has considered the initial state [176162117711781179118013911811182] 



ip{x,t = 0) = (e' 



— ik2. 



)Q{-x), 



(81) 



which continuously vanishes at a; = 0. It follows then, that the time-dependent solutions 
along the internal and transmitted regions of the potential become, respectively. 



ipix, t) = (/)(x, k)M (0, fc, t) - (t>{x, ~k)M (0, -fc, t) 
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Fig. 17. Comparison of the wavefronts of the free-wave evolving solution \il>o{x, t)p and the transmitted 
solution at resonance energy \tp{x,t)\'^ through a double-barrier system at a long time t= 1540 ps, to 
show that the delay time is of the order of 2h/ri , with Fi , the corresponding resonance width. 

oo 

- Mx)M{0,kn,t); {0<x<L), (82) 

n— — QC 

and 



ipix, t) = T{k)M (x, k, t) - T{-k)M (x, -fc, t) 

OO 

- TnM{x,kn,t)- {x>L), (83) 

n— — oo 

where the terms (f>n{x) and Tn are now given by (pn{x) — 2ikun{0)un{x) / {k'^ — fc^) 
and Tn = 2iA:M„(0)u„(L)e~*'^"^/(fc^ — fc^J. It may be shown that at very long times 
M{x,—k,t) goes to zero [61] and hence the above solutions tend also to the stationary 
values given by Eqs. (HH) and ([50)1 . 

In [391 an alternative procedure was derived to obtain the time-dependent transmitted 
solution. This procedure starts from the expression for the time-evolved wave function 
along the transmitted region 

oc 

^(l){k')T{k')e''''''-''"'"'/'^"', (84) 

— oo 

where 0(fc') is the A:'-space wave function (i.e., the Fourier transform of the initial wave 
function) defined by 

oo 

^{k')= J -^e-'='-V'(x,0). (85) 
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T(fc) for fc < is the analytical continuation in that domain of the transmission amplitude 
for left incidence, T{k > 0), not to be confused with the transmission amplitude for right 
incidence 60]. 

For the initial state given by Eq. (|8ip. it follows from Eq. ([85|) that 



oo 




_ i / 1 1 
~ VSvr \k' - k + ie k' + k + ie 

where e is an infinitesimal positive number. The other ingredient to evaluate Eq. (j84|) is 
to make use of the Cauchy expansion for the transmission amplitude derived in Appendix 
m i.e., Eq. (|K26)) . written as 

n— — oo ^ ^ 

where r„ = iu„(0)u„(L) exp(— ifc„L). Substitution of Eqs. ([86|) and ([87|) into Eq. (|84l) 
leads, after some algebraic manipulation, to exactly the same expression given by Eq. 
(|83p . This procedure allows one to obtain more easily the time-dependent solution than 
the Laplace transform approach since it facilitates the consideration of other initial states. 
We shall refer to this approach as the Fourier transform approach. Its extension to deal 
with the internal region is straightforward [184j . 

There have been two main lines of research where the above exact time-dependent so- 
lutions for the wave function have been used in the investigation of transient phenomena. 
One line has addressed detailed studies on the dynamics of resonant tunneling on and off 
resonance energy in double |61ll76ll86ll87j . triple [179] and multiple |177ll88j resonant 
tunneling structures, whereas the other line of research, considered in the next section, 
has addressed the issue of transient effects in relation to the tunneling time problem 
|()2ll78l()3ll8()i:fflll74ll8lTT82] . 




6.1.1. Buildup dynamics of the transmission resonance spectra 

Most studies involving the quantum shutter setup consider a given value for the energy 
of the initial state, and evaluate the time-dependent solution for the probability density 
either as a function of time for a fixed value of the distance or as a function 
of the distance for a fixed value of the time. In Ref. [177j . however, Romo addresses the 
issue of the building up of the transmission resonances in a superlattice of finite range 
by considering the time-dependent solution for fixed values of both the distance and the 
time and varying instead the energy of the initial state in a relevant energy range. He 
considered the initial state given by Eq. ([5T|l which leads to the time-dependent solution 
given by Eq. (j82p to study this transient behavior for different superlattice potential 
profiles. In a finite superlattice involving -I- 1 barriers the resonances group themselves 
in minibands involving N resonance levels isolated from each other. Romo studied the 
building up dynamics for a given miniband and found that the solution for |';/'(a;,t)p 
using Eq. (j82p may be written as a simple analytic expression which involves only the 
N resonance levels of the miniband. Figure [18] exemplifies his findings. It consists of a 
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Fig. 18. Snapshots of the time-dependent probability density for a periodic superlattice of five barriers, 
with parameters given in the text, reproducing Fig. 4 of Ref. |177| . The values of \'^(E, L; t)p vs E are 
calculated Eq. (7) of that paper at different fixed times (solid lines); (a) t=0.6 ps; (b) t=2.0 ps; (c) t=5.0 
ps; (d) t=80.0 ps. The transmission coefficient is also included for comparison (dotted lines). 

superlattice formed by five rectangular barriers with heights Vq — 0.2 eV and widths 
feo = 5.0 nm and well widths Wq = 5.0 nm, and it shows a series of plots of | 'i/^{E, L; t) | 
as a function of the energy E at different fixed times. In the snapshots depicted in 
Figures [TSTa) to [TSl d) one appreciates the "birth" of the transmission resonances and 
their subsequent evolution towards the stationary regime. At the very beginning of the 
tunneling process, one appreciates a smooth curve with no peaks, see Fig. llSf a): that 
is, no evidence of the resonances have appeared at this early stage. However, as time 
elapses some peaks in the curve gradually begin to appear as can be appreciated in Figs. 
[rST b'l to [TSTd) . The height of the resonance peaks increases at different rates towards 
their corresponding asymptotic values. Notice that the buildup of the peaks is faster 
for the wider resonances {n = 2 and 3), and slower for the thinner ones {n = 1 and 
4). It is also obtained that the buildup of the transmission peaks is governed by the 
analytic expression T'P'^°''[1 — exp(— t/i^)], where T^'^"'^ is the height of the corresponding 
transmission peak and = 2fi,/r„, namely twice the lifetime fi/r„ of the resonance level. 



6.2. Wavepacket initial states 

The transmission of wavepackets through one-dimensional potentials is a model that 
has been of great relevance both from a pedagogical point of view, as discussed in many 
quantum mechanics textbooks, and in research, particularly since the advent of artificial 
semiconductor quantum structures |207l208j . In [182] . Yamada, Garcfa-Calderon and 
Villavicencio extended the quantum shutter initial condition to certain type of wave 
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packets. They considered the foUowing initial condition, 



^(^'0)=^ / '^M fc_fc„^,^ +c.c.), (88) 



where A = -^/A {1 + (A/fco)^}/27r with A > 0, and c.c. stands for complex conjugate. An 
important feature of this initial state is that it automatically vanishes for x > 0. This is 
immediately seen from the fact that the integrand e*'^^/ (fc — /cq + i^), which corresponds 
to a Lorentzian momentum distribution centered at hko with width HA, has a simple 
pole only in the lower-half of the complex fc-plane. An explicit expression for Tp{x, 0) can 
be easily obtained by the method of residues. It is found that 

A-kAc^"" sinfeox, x < 
■4){x]t = 0) = { . (89) 

0, a;>0, 

A measure of the packet width is 1/A. It can be easily proved that the wave packet is 
normalized, i.e., J da;|?/;(a;, 0)p = 1. The Fourier transform of the above expression is 

m = ^/2^A ( ^ ^—-\ , (90) 

and then, following the same procedure of the previous example, these authors arrive at 
the solution 



4,{x, t) = -^ VA{1 + (A/A:o)2} 



r(fco - iA)M{x,ko - iA;t) 



- T(-fco - iA)M {x, -fco - iA; t) 



2ko J2 



2 

n— — oo 



fco - (kn + «A)- 



:M{x,kn;t) 



(91) 



where, we recall that r„ — iu„(0)it„(i) exp(— iA:„L). 

Most time-dependent numerical studies consider Gaussian wavepackets as initial states 
|208l209l210l211j . though in some recent work, the formation of a quasistationary state in 
the scattering of wavepackets on finite one-dimensional periodic structures also involves 
also some analytical considerations [212]. As discussed above, analytical approaches have 
been mainly concerned with quasi-monochromatic initial states in a quantum shutter 
setup. In some recent work, however, analytical solutions to the time-dependent wave 
function have been discussed using initial Gaussian wavepackets for square barriers [200j , 
delta potentials |213j and resonant tunneling systems near a single resonance |214j . 

In Ref. [185], Villavicencio, Romo and Cruz derive an analytical solution to the time- 
dependent wavefunction for an initial cutoff Gaussian wavepacket, for the free evolving 
case and two barrier potentials: the (5-potential and a very thin and very high rectangular 
barrier. They consider the Fourier transform approach and the analytical solution follows 
provided the center of the Gaussian wavepacket is far from the interaction region, as is 
usually assumed on physical grounds. In such a case, the tail of the Gaussian is small 
near the interaction region, a regime they refer to as the small truncation regime. 
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The cutoff Gaussian wave packet is given by the expression 



M^) = { , (92) 

0, x>0 

where Aq is the normaHzation constant, and xq and a stand, respectively, for the center 
and effective width of the wavepacket. The interesting point of this approach is that 
substitution of Eq. ([92| into Eq. ((85|l gives an exact expression for the corresponding 
Fourier transform, 

Mk) = Aaw{iz), (93) 

where Aq = {l/V2^)[i2TTa'^)^/'^/y/u{izo)], z = a:o/(2a) - i(k - fco)cr, zq = xo/(\/2(j), 
and u'(y) is the w-function [75] (see Appendix [C|) . Along < and z 3> 1 using the 
relationship w(iz) = 2e^ —l/(7r-'^/^z) — l/(27r-'^/'^z'^) + ..., allows to write w(iz) ~ 2exp(z^). 
Substituting this last expression into Eq. (|93p and the resulting expression into Eq. ([84|l 
yields 

/dV 2 / ,2 

v 27r 



For T = 1, the above expression leads to an analytical expression for the free evolving 
Gaussian wavepacket |185j . 

e,p(- ["-;°-(^V^)^]n , (95) 



40-2 [i + it/r] 

where r = 2ma^/h, that is identical to the analytical solution for the extended free 
Gaussian wavepacket, i.e., no cutoff Gaussian wavepacket. The above authors considered 
the quantum shutter setup with a cutoff Gaussian wavepacket initial state for a delta- 
barrier potential of intensity A, i.e., V{x) = XS{x). The corresponding transmission 
amplitude reads 

Notice that the transmission amplitude given above has just one pole, an antibound pole 
at ka = —i{m\/h^). Substitution of Eq. ((85|) into Eq. ((84|) leads to an expression that 
acquires a simple analytical form in the small truncation regime, namely, 

i^{x, t) = i^i{x, t) - C7e*fco^-^?ifco2t/2m^(^/^ ^/^ (97) 

where '4'i{x,t) stands for the free solution, C = (7r/2)i/'*(2mcri/2A)/n2 ^nd M{x',q',t) 
stands for a Moshinsky function where x' = x — xq — vot, q' = — (fco + im,\/h?), and 
t' = t — ir. The solution given by Eq. ((97|) holds also for a rectangular barrier of height 
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(a) 




Fig. 19. Figures (a), (b) and (c) display (T|i/'(a;, t)p as function of time in units of the lifetime = fi/^i 
(full lines) at different distances from the DB system, calculated using Eq. H99| l: (a) 2L, (b) 200L and 
(c) 2 X lO^L, where L is the length of the system. Also shown in these figures is the free evolving cutoff 
Gaussian wavepacket (dotted line). Figure (d) exhibits the transmission spectra of the DB system as a 
function of energy, in units of the barrier height, using the exact numerical calculation (full line) and 
that obtained using the resonance expansion given by Eq. l|E.26| l (dot line). See text. 

vq and width L in the hmit of small opacity [2mVa]^/^ /h <C 1, by letting A = VoL and 
performing the translation x — > {x — L). 

Recently, Cordero and Garcia-Calderon [215j have obtained a formal solution of the 
transmitted time-dependent function for an arbitrary potential profile as an expansion 
in terms of resonance states and poles of the problem that follows by substitution of Eq. 
(|R26| into Eq. (gll), namely, 



1 



OO 

E 



-iknL 



dk 



(l)o{k)e 



ikx — ihk i/2m 



(98) 



where p„ = u„(0)u„(L)/fc„. For a cutoff Gaussian initial wavepacket, substitution of 
Eq. into Eq. leads, in the small truncation regime, to the expression 



V(x,0=T(fco)V'f(x,t) 
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1 

kn ^0 



(99) 



where D{zo) = — 2i(27r)^/''-\/o'/erfc(zo) and the Moshinsky function M(x' ,t';k'^), de- 
pends now on fc„ — kn — ko, with fc„ the complex resonance pole. The other parameters 
remain as defined above. Figures [19] (a)- (c) exhibit a\ip{x,t)\'^ in units of the lifetime, for 
a cutoff Gaussian initial wavepacket with parameters xq = —5.0 nm, a = 0.5 nm and an 
incident energy Eq of half the potential barrier height V, tunneling through a double- 
barrier resonant system (DB) with the same parameters as given in the discussion in 
16.11 The profile of the transmitted wavepacket exhibits a transient behavior that after a 
very large distance, acquires a fixed shape that resembles, in time domain, the resonance 
spectra of the system. As a comparison, Fig. [19] (d) yields the transmission coefficient 
of the DB system vs energy in units of the potential height, using an exact numerical 
calculation (full line) and the resonance expansion given by Eq. ([KM]) . 



6.3. Relativistic effects in quantum transients 

An interesting feature of the solutions for cutoff initial waves, occurring both in the free 
case |183j and in the presence of a potential interaction [HT], is that, if initially there is a 
zero probability for the particle to be at x > 0, as soon as i 7^ 0, there is instantaneously, 
a finite, though very small, probability to find the particle at any point x > 0. In Ref. 
|31j . Garcfa-Calderon, Rubio and Villavicencio, studied the short time behavior of the 
solution of the quantum shutter in the presence of a finite range potential extending 
through the region < a: < L with the initial conditions given, respectively, by Eqs. (|76p 
and (|8T|) . which lead along the external region, x > L to the solutions given by Eqs. (jTS)) 
and (|83p . The found that at short times the above solutions behave as 

ip{x,t) ~ -t^/^; x>L. (100) 

X 

where A is a constant. The above expression tells that the probability density at any 
distance x from the potential will rise instantaneously with time. Although this should 
not pose any conceptual difficulties because the treatment is not relativistic, the above 
authors ask themselves whether the above nonlocal behavior arises because the initial 
condition is a cutoff wave. In order to solve this question the authors considered the 
solution of the Klein-Gordon equation for a delta potential V{x) = bsS{x) by using 
Laplace transforming techniques and found an analytical solution of the problem that 
at short times becomes different from zero after a time to = L/c, with c the velocity of 
light, thus restoring Einstein causality. Hence the non-local behavior of the Schrodinger 
solution is due to its nonrelativistic nature and not a result of the quantum shutter setup. 
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7. Time scale of forerunners in quantum tunneling 

Quantum tunneling, the possibility that a particle traverses a classically forbidden re- 
gion, constitutes one of the paradigms of quantum mechanics. Most textbooks usually 
address its stationary aspects only, and solve the stationary Schrodinger equation at a 
fixed energy E. In the time domain the analysis of the transient behavior of the time 
evolution of the wave function along the internal and transmitted regions of a potential 
has been considered in the framework of the quantum shutter setup as a natural exten- 
sion of the free case considered by Moshinsky. Since the seminal work of Biittiker and 
Landauer |189| . the time-dependent aspects of tunnehng have attracted a great deal of 
attention. Much of this work has led to controversy, as several authors have proposed 
and defended different "tunneling times" . In fact each has its own virtues, weaknesses, 
physical content, and range of applicability. For reviews see |190|191|T92 193 194j. 

An analysis of the involved non-commuting observables (the projectors that determine 
the final transmission and the probability to find the particle in the barrier region) shows 
that, from a fundamental perspective, there is no unique tunneling time, because several 
quantizations are possible due to different operator orderings and defining criteria [195j . 
It is thus necessary to specify precisely how to time the quantum particle in the tunneling 
regime, since different procedures lead to different relevant time scales. For example, the 
traversal time of Biittiker and Landauer (BL time) [189j . tbl ^ L/vsc, given by the 
barrier length L divided by the "semiclassical" velocity Vsc — [2(V^ — E)/m]^/^, marks 
the transition from sudden to adiabatic regimes for an oscillating barrier |189j . and 
determines the rotation of the spin in a weak magnetic field in opaque conditions [196j : 
whereas the average over wavepacket components of the (monochromatic) "phase times" 
provides the mean arrival time of the transmitted wave packet [1951197] . These two 
time scales may be very different. The Biittiker-Landauer time increases with decreasing 
energies up to a finite value, whereas r^''(0, L) (the so called extrapolated phase time) 
diverges as £' ^ 0, and tends for increasing i to a constant value, 2h/[veci'2mEy/^]. This 
latter property implies that the arrival of the transmitted wave becomes independent of 
L (Hartman effect |198j ). although the independence only holds until a certain critical 
length Lc |195j where above-the-barrier components start to dominate. For L > Lc 
the mean arrival time depends on L linearly. While tbl and r^'* are surely the most 
frequently found tunneling times, they do not exhaust all timing questions as we shall 
see. 

7.1. Transients for a square barrier 

A simple and physically interesting complication of the elementary shutter is the addi- 
tion of a square barrier after the shutter edge. Brouard and Muga studied this problem 
paying attention to the "tunneling" configuration in which the carrier energy is below 
the barrier top, Eq < Vq [60J. They used a complex momentum plane approach asso- 
ciated with asymptotic expansions which has been applied later to different transient 
phenomena |199l63l64l32l200l36l65l201l202j . The aim of the approach is to describe the 
wave propagation with a minimum of elements and maximal efficiency. It provides the 
wave function, as well as characteristic velocities and times, and localizes the origin of 
the main contributions at critical points in the complex energy or momentum planes 
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such as poles, saddle points, and branch points, following the analysis by Sommerfeld 
and Brillouin for the propagation of light in dispersive media j203| . This is useful con- 
ceptually, and improves the efficiency of the numerical calculations. Pioneering work in 
this direction was done by Stevens [58^, who proposed a sequence of "tunneling prob- 
lems" for the step potential barrier. Using approximate arguments, his main conclusion 
was that for energies below the barrier height Eq < Vq, a, wave front traveling with the 
semiclassical velocity Vsc = [2(Vb — Eo)/m]^^'^ could be identified. He related this wave 
front to the crossing of a pole by the steepest-descent path. This conclusion has been 
later shown to be generally unjustified, due to the simultaneous effect of other critical 
points |171|59|115|116j . More on this in l7. 21 below. The study of finite width barriers is a 
natural extension of Stevens's tunneling problems for the step potential. Along this line, 
Jauho and Jonson 5S] examined numerically the propagation of an initially sharp packet 
(a cutoff plane wave) through a square barrier, and Moretti carried out an approximate 
asymptotic analysis for this system IIX well. 

Brouard and Muga [60 provided an exact solution of the dynamics that retains the 
basic philosophy of working out the time-dependent wave function by contour deforma- 
tion in the complex momentum plane and extracting contributions from critical points, 
which in most cases reduce to a w-function. 

It is instructive to note the main features of the analysis in [60] by first studying 
Mosliinsky's simplest shutter problem. The general solution for free motion takes the 
form 



Mx,t)= J {x\k)e-'^'''/''{k\ip{t = 0))dk, (101) 

— CSO 

where Ek = k^h'^/2m and {k\k^) — S{k — k'). For an initial state corresponding to a 
cut-off plane wave e*^"^B(— x), 

OO 2 / s 

^ r ^i{kx — k ht/2m) 

— OO 

The steepest descent path (SDP) is a —45° straight line cutting the real k axis at the 
saddle point (fc ~ mx/th) of the exponent. For a fixed position x > 0, the saddle point 
moves from oo to as i grows from to oo. Deforming the contour along that path, 
completing the square, and using the same u- variable of Eq. ([5]), 

fcj-\ 1/2 

nt \ / , mx \ 



which is zero at the saddle, and real along the SDP, we get 

^imx^ /2th r g"'^^ 

V'o(a;,t)=«- / du {x,t>Q), (104) 

27r J u — Uq 

where mq ~ u{k = fco) is the pole position in the u-plane. Note that it "moves" with t 
and/or x; by contrast, the pole in the fc-plane is fixed at fco- The contour r„ goes along 
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the real u axis from — oo to oo plus a circle around uq whenever its imaginary part is 
positive. Finally, 

g'ima;^ /2th 

Mx,t)^ wi-uo) ix,t>0), (105) 

which is Moshinsky's solution M{x,k,t), Eq. (|4]). At fixed x, Uq ^ ik^ht/2m as t ^ oo 
and a stationary regime with constant density is reached. At small times the saddle at 
w = is far from the pole and it is the only important critical point. Setting u = in 
the integrand, 

W-'OP-^;^ (106) 

The treatment for the square barrier starts from an integral like (jlOip . substituting the 
plane waves by scattering states. Thus, for the transmitted part, x > a, 

oo 
— oo 

which is in fact a generic result for cut-off potentials and initial states. The difference 
with the elementary treatment of the free motion case is the need to consider the complex 
poles of T(k) in the lower half-Zc-plane, kj, which leads to a series in terms of w-functions 
for each of them. The short time asymptotic behavior does not change with respect to the 
free case, Eq. (|106p , because the poles are yet very far from the saddle and integration 
contour. As time progresses the saddle gets near the poles, and a velocity of propagation 
can be assigned to each resonant pole contribution from the condition of the crossing 
of each pole by the contour. These contributions are transient, fade away with time, 
their residue terms being proportional to g~'^k^ht/2m ^ only the contribution of the 
structural pole fco remains eventually^ which sets in around the phase time. 

The internal part, < x < L, becomes more complicated in this treatment because of 
the interference between different terms, but the results show clearly that a front moving 
with a semiclassical velocity is not seen in the exact solution. 

Alternatively, using the exact analytical approach discussed in the previous section, 
Garci'a-Calderon [206] examined the time evolution of the initially cutoff wave given 
by Eq. ([5T|l with energy E — 0.01 eV on a rectangular barrier. The parameters of the 
barrier are typical of semiconductor quantum structures [207) : Barrier height Vb = 0.23 
eV, barrier width L = 5.0 nm and effective electronic mass m = 0.67me, with me the 
bare electron mass. The results of the calculation exhibit an absence of a propagating 
wave along the internal region of the potential. The absence of a propagating wave along 
the tunneling region may also be inferred from the corresponding analytical solution (see 
Eq. (|82p since the Moshinsky functions there arc independent of x. The inset to Fig. [201 
exhibits the probability density as a function of time at a fixed internal distance xq — 2.0 
nm. A peaked structure at t ~ 5.6 fs is clearly appreciated. This is a replica of a similar 
structure that appears at a; = L, named time domain resonance 62 178 180 which is 
further discussed in l7.3l and l7.4l Similar non-propagating structures with different peaked 

* The "structural poles" are associated with the structure of the wave-function and are to be distin- 
guished from poles of the resolvent, in particular from resonance poles |204l205l . 
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Fig. 20. Plot of the probability density normalized to the incident flux along the internal region of a rect- 
angular potential, denoted by \tf'(x,t)\^, for various times which exhibits that there is no a propagating 
wave along the internal region: ti = 1.0 fs, t2 = 3.0 fs, and = 30.0/s. For comparison the stationary 
solution has been also included (dashed line). The inset shows a replica time domain resonance at the 
internal distance xq = 2.0 nm. See text. 
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Fig. 21. Exact calculation of tp (solid dots) as a function of the position x. The full square indicates the 
value of the penetration length . See text. 

values tp appear at every internal point. In Ref. [63] this was further investigated. Figure 
\n] exhibits a plot of tp (solid dots) as a function of position along internal values of 
the distance xg of a rectangular barrier with parameters: Vq — 1.0 eV and L = 40.0 
nm. Figure TIT\ explores distances of the order of, or smaller than, the penetration length 
defined as Kq^, where kq = [2rn(Vb — Eo)]^^^/h, and Eq — fi^A:Q/2m stands for the 
incident energy. 

One might clearly identify two regimes. The first of them corresponds to a basin, 
where for small values of the position, tp decreases reaching a minimum value as x 
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increases. However, if x is further increased, tp begins to grow as a function of the 
position. Apparently this second regime corresponds to a situation where tp increases 
Unearly with x. The resuh depicted in Fig. [21] appears also to be a replica of a similar 
behavior of the probability density at the barrier edge x = L, reported in Ref. [62] ■ Before 
we further clarify and interpret these effects in 17.31 it will be useful to step back from 
the relative complexity of the square barrier and analyze an even simpler case. 



7.2. Muga-Biittiker analysis of evanescent waves 



Here we shall follow [M] to study the role of the Biittiker-Landauer time in quantum 
transients due to the injection of particles with energy below the potential level, so that 
an evanescent wave is finally formed. For the Schrodinger equation 

with V ^constant, we define wavenumber k and frequency lo as 

uj = Elh=V/n + k^h/2m, (109) 
k=[2m{uj-V/K)/hYl^. (110) 

Particle injection is modeled by fixing the wavefunction at the origin for all times, 

'ip{x^{),t) = Q{t)e'"^°\ (111) 

for a frequency loq — hkQ/2m, and assuming ip{x,t < 0) = 0. 

The inverse Fourier transform of the source "00 (2; = 0, /coi — e~"^"*0(t) is given by 

00 

^{x = 0,ko,Lo)^^ [ M^-O,ko,t)e'^'dt^^ ^— . (112) 

V27r J V27r cj - cjq + ^0 

— oo 

For X, i > 0, the wave function evolves according to 

oo 

^ r ^ikx—iujt 

Mx,ko,t) ^ — / du; — , 

— oo 

where Imfc > 0. If we are interested in carrier frequencies below threshold, loq < V/h, 
it is also useful to define kq — [2m{V/h — ujq)/K\^/'^. The integral can be written in the 
complex fc-plane by deforming the contour of integration to r_|_, which goes from — oo to 
oo passing above the poles, 

j /• ^ikx-i{k^h/2m+V/h)t 



ipoix,ko,t) = — I dk2k 



r 



+ 



* / I \ ^ikx-ik-'ht/2m-iVt/h 



2tt J \k + ino k — iKo 
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The contour can be deformed further along the steepest descent path from the saddle at 
ks = xm/th, the straight line 



xm/th (113) 



(k ^ kji + ikj), plus a small circle around the pole at fcp = iKQ after it has been crossed 
by the steepest descent path, for fixed x, at the critical time 

Tc = xfn/h[Re{ko) + Im(fco)]. (114) 

For tunneling Re(fco) = so that Tc — tbl, where now 
xm 

is a slight generalization of the Biittiker-Landauer time, with x playing the role of L. 
The solution, using the -u- variable (I103p . real along the steepest descent path, is 

^ — iVt/h^ix^'m/2th 

^ix,t) = [w{-uo) + wi^u'o)] , (116) 

where uq = u{k = /cq), and Uq = u{k = — /co)IZ] Eq. (|116p can be written for xkq 3> 1 
(opaque conditions) as the sum of contributions from the saddle point and the pole, 

ij(x,t)^^p{x,t) + i;s{x,t), (117) 
i;p{x, t) = e— ot+^fco.Q _ ^ (118) 

M^,t)^ — — _ + _ . (119) 

ZITT'-'^ \Uq UqJ 

Asymptotically, only the "monochromatic" pole contribution remains and a stationary 
evanescent wave is formed, but, before that regime, a forerunner (bump in the density 
with a well defined maximum) is observed at fixed x. For carrier frequencies below the 
cut-off tb l plays clearly a role in the exact solution since 

1 + i ht , , , 

uo^——\—HQ(i-TBLt), (120) 
2 V ni 

, 1 + i ht 

Uo = ^-^ —Koi-l-TBL/t). (121) 

However, it cannot be identified with the transition to the asymptotic regime since the 
saddle contribution dominates up to an exponentially large time ttr = Ur (x) > tbl that 
can be identified from the condition \tpp\ = \ips\ [Si]- (By contrast above threshold, Tc 
in Eq. (|114p is a good scale for the arrival of the main peak.) A different instant is the 
time of arrival, tp = tp{x), of the temporal peak of the forerunner at a point x. It can 



^ This formal result holds for all carrier frequencies uJo, above or below threshold. Putting V = it reads 
simply jI>{x, t) = M{x, ko,t) + M{x, — fco, t) in terms of Moshinsky functions, i.e., the source boundary 
condition imposed at a; = gives the same wavefunction as the shutter with reflection amplitude R = 1 
opened at t = 0. The case R = —I and its relation to inhomogeneous sources was examined by Moshinsky 
[216] and is reviewed in Appendix [P] 
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be obtained from from d\%lJs\'^ /dt = and, for opaque conditions, tp turned out to be 
(surprisingly) proportional to the BL time tbl , 

tp = TBL/i^'^, (122) 

even though the time-frequency analysis of the forerunner |217l34j confirmed that it was 
composed by frequencies above threshold so it was not tunneling[£] The frequency of 
the saddle is ujs ~ Eg/h, where Eg ^ V + x'^m/{2t^) ^ V + k^fi^ /{2m), coincides with 
the energy of a classical particle traveling from the source to a; in a time t. It is very 
remarkable that in spite of its frequency content, the forerunner's peak "travels" with a 
velocity proportional to Vsc, Vp — S^^^Vsc, which increases with decreasing energies, and 
its intensity diminishes exponentially as it progresses along the coordinate x. Other works 
had already pointed out the dominance of non-tunneling components in the forerunner 
[218|219|171|59|60] but had not characterized its time dependence. 



7.3. Forerunners which do tunnel 



Biittiker and Thomas |33| proposed to enhance the importance of the monochromatic 
front associated with ujq compared to the forerunners by limiting the frequency band of 
the source or of the detector; it was shown later in [34] that the monochromatic front 
could not be seen in opaque conditions even with the frequency band limitation. How- 
ever, a clear separation of the amplitude into two terms, one associated with saddle and 
forerunner and the other with the pole and the "monochromatic front" , is only possi- 
ble for opaque conditions. Non-semiclassical conditions have been much less investigated 
[220j . even though these are actually easier to observe because a stronger signal may be 
obtained. Ref. [63j and this subsection are mainly concerned with them. 

Garcfa-Calderon and Villaviccncio [62] examined the time evolution of an initial cutoff 
wave truncated at the left edge of a rectangular barrier potential (quantum shutter setup). 
There it was found that the probability density at the barrier edge x = L, exhibits at 
short times a transient structure named time domain resonance, see Fig. 1221 The behavior 
of the maximum, tp, of the time domain resonance showed a plateau region for small L, 
or more accurately a shallow basin, followed by a linear dependence for larger L; this 
behavior is reminiscent of the Hartman effect, but the time of the plateau did not coincide 
with the phase-time estimate. It was found that for a broad range of parameters, tp in 
the basin may be written approximately as 

*p = 7^' (123) 

where ei and Eq, correspond to the energy of the first top-barrier resonance and the 
incidence energy, respectively. On the other hand, along the linear regime, at larger 
values of L, tp is described by 

t^p = ^, (124) 



^ The simplest time-frequency information is contained in the local instantaneous frequency. Writing 
ip{x,t) = |'!/>(x, t)|e"^(^''' it is defined as ij{x,t) = —d<f>{x,t)/dt. More sophisticated treatments such 
as spectrograms, Wigner functions or other time-frequency distributions provide further moments and 
information |217l34j . 
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Fig. 22. Time evolution of \jp{x,t)\'^ for a rectangular barrier of height V = 0.3 eV and width L = 5.0 
nm, evaluated at a; = L, exhibiting a time domain resonance peaked at tp = 5.3 fs. The incident energy 
is E = 0.01 eV and the effective electron mass is m = 0.067me, where rrie stands for the electron mass. 
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Fig. 23. Maximum of the time-domain resonance tp (solid dots) as a function of the barrier width L with 
fixed barrier height V = 0.3 eV and incidence energy E = 0.001 eV. For comparison, are also plotted 
the Biittiker traversal time ts, the Biittiker-Landauer time tsL and the phase-delay time to. 

where vi — hai/m, with ai the real part of the first top-barrier pole ki = ai — ibi. The 
above considerations are illustrated in Fig. 1231 This figure also shows a comparison with 
exact calculations for the Biittiker traversal time t b , the semiclassical Biittiker-Landauer 
time tsL and the phase-delay time t^, see Eqs. (3.12), (1.7) and (3.2) of [196]. Notice 
that Ib and Ibl do not describe the basin regime and that they are close to the values 
of tp along the regime that increases linearly with L. 

It was shown in [63j that this basin dependence, and the corresponding time scale, is 
not only present at the transmission edge of the square barrier studied in [62J . It may also 
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be found, mutatis mutandis, in the sharp onset source model with constant V examined 
in for small xkq. As discussed above, the basin is also found for the time of the 
forerunner versus position in the internal region of the rectangular barrier, and for a step 
potential barrier. There are no resonances in the sharp onset source model, or for the 
step potential, so the time scale of the forerunner at the basin minimum is in these cases 
inversely proportional to Kg, namely to the difference between the "potential level" V 
and the source main energy Eq. This is to be contrasted with the dependence on Kq ^ of 
the traversal time tbl- Apart from certain peculiarities, all models show a small length 
region of the order of the penetration length where the forerunner is dominated by 
tunneling components and arrives at a time proportional to Hq"^ ■ 

The time-frequency analysis showed that the peak of the forerunner at small lengths, 
in the evanescent region, is composed predominantly by under-the-barrier components, 
so that indeed there is a genuine tunneling time scale different from the phase or BL 
times [63] , 

In Ref. |178j the analytic solution to the Schrodinger equation for cutoff wave initial 
conditions given by Eq. (|83p was used to investigate the time evolution of the transmitted 
probability density for a rectangular potential. For a broad range of values of the opacity 
defined as 

a = [2mV]^/^L/h, (125) 

it was found that the transmitted probability exhibits two evolving structures. One 
refers to a forerunner which follows from the propagation of the time-domain resonance 
whereas the other structure consists of a semiclassical propagating wavefront, that re- 
sembles the time evolution of the free term modulated by the value of the transmission 
coefficient. A regime is also found for values of the opacity a below a critical value ac 
where there are no time domain resonances and hence no forerunners. Moreover, it is 
also found that this situation is related to a time advance, i.e. a positive delay time. This 
occurs for very thin or very shallow potential barriers. For the effective electron mass 
m = 0.06777ie, typical of GaAs semiconductor materials, ac = 2.0653. 

7.4. Barrier opacity and time scale for tunneling as a transient effect 

The short-length tunneling time scale discussed before is further investigated in Ref. 
jl80j : Garci'a-Calderon and Villavicencio look at the peak value of the probability density 
at the barrier edge a; = L in the shutter problem and characterize the regime that 
corresponds to under-the-barrier tunneling by the opacity of the system, i.e., Eq. (|125p . 
A time-frequency analysis establishes the existence of a range of values of the opacity, 
independent of the incidence energy, where the peak value of the time domain resonance 
is governed by a single frequency, which suggests that the system acts as a frequency 
filter. Figure [24l exhibits the relative frequency tOav/^v as a function of the opacity for 
three different values of the parameter u ~ V/E. Here ujav = Im[((i'0/di)/V'], and the 
cutoff frequency LOy = V/H. The time-dependent solution is given by Eq. ([5^ . Notice 
the existence of an opacity "window", in the range of values 2.065 < a < 3.3, where the 
average frequencies are always below the cutoff frequency, i.e.. In this example the value 
oi V = 0.3 eV and the mass m — 0.067r7ie, with rUe the bare electron mass. 
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Fig. 24. Relative frequency LOav/i^v measured at the barrier edge function of the opacity 

a. The parameters: u = 300 (sohd dot), u = 10 (sohd triangle), and u = 5 (solid square). Note that 
for values of the opacity smaller than a ~ 3.3, the relative frequencies for all values of u are below the 
cutoff-frequency uJav/i^v = 1 (dashed line). See text. 

7.5. Time scales for relativistic equations 



In [32 the authors investigated the role played by tbl, and the characteristics of 
the forerunners at large and small distance from the source for under-cut-off, unit-step 
modulated excitations with the (dimensionless) relativistic wave equation 



■ 



* = 0, (126) 



which turns out to be more accessible experimentally (with waveguides in evanescent 
conditions |220|221j ) than the corresponding Schrodinger equation 

id'^/dt= -d'^^/dY^ + (127) 

The stationary equations corresponding to Eqs. (|126p and (|127p are equal and have 
equal solutions, but the two cases are not equivalent in the time domain, because of the 
dispersion relations between frequency and wavenumber. In fact a wealth of qualitative 
changes with respect to the Schrodinger scenario were found [32] , 

Eq. p26|) leads to the Klein-Gordon equation for spin-0 particles by means of the 
substitutions 

Y = x{moc)/h, (128) 
T = tmoC^/h, (129) 

where x and t denote dimensional position and time, and mg and c are the particle's rest 
mass and the velocity of light in vacuum. Another interesting connection, much more 
promising to implement the analysis experimentally [221] . and free from the conceptual 
puzzles of the former, may be established with the equation that governs the electro- 
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magnetic field components in waveguides of general constant cross section with perfectly 
conducting walls. This requires the substitutions 



Y = x\, 
T = t\c, 



(130) 
(131) 



where A > is one of the eigenvalues of the waveguide. The main differences between 
the relativistic and non-relativistic cases are: (a) the relativistic solutions are not simply 
related to each other by time and position scaling as in the non-relativistic case, so that 
qualitatively different "shallow" and "deep" tunneling regimes may be distinguished; (b) 
tunneling is more robust relativistically, both in the precursors (more tunneling frequen- 
cies) and asymptotically (more "density"); (c) the "first" relativistic precursor, right 
after the limit imposed by causality does not have a non-relativistic counterpart, and 
does tunnel; (d) the "second" precursor, which tends to the non-relativistic one for ex- 
citation frequencies near cut-off, has an oscillating structure that may be characterized 
by its envelopes. The traversal time tbl is not an exact measure of their arrival except 
in the non-relativistic limit. 

8. Transient phenomena caused by a sudden perturbation 

We now turn our attention to the dynamics following a general sudden perturbation 
in which the potential of the system changes from V- (x) to V+ {x) abruptly. 

8.1. Potential steps, resonant tunneling 

One of the most important tunneling processes is resonant tunneling, which has been 
employed in several semiconductor devices since the original double-barrier diode of Tsu 
and Esaki |222j . and more recently to control the energy levels, transport through, and 
state occupancy of semiconductor quantum dots and wells [223| and nanowirc resonant 
tunneling diodes |224j . Quantum dots are candidates as elements in quantum circuits, and 
basic quantum manipulation is possible by applying sharp voltage pulses or microwave 
signals to the contact leads surrounding the dot region 225.226,227^ . The response of 
quantum dots to these stimuli in the sub-nano-second to almost milli-second time domain 
is of much current interest. Two basic manipulations that can significantly affect charge 
transport across the quantum well (dot) are the change (i) va. the source-drain voltage 
applied across the quantum well (dot), which tilts the whole potential structure, and 
(ii) in the gate electrode voltage, which may modify essentially (and electrostatically) 
the quantum well (dot) depth [228 229j. At resonance, full transmission occurs, whereas 
ofF-resonance, there is essentially full reflection. Understanding the transitions that fill 
or deplete a resonance or a bound state following a change in potential and determining 
the speed of the change is clearly of fundamental interest. 

For scattering systems, an elementary transient process is the transition between eigen- 
states of the original and final Hamiltonian in the continuum. A stationary scattering 
wave function responds to an abrupt change in the potential shape by forming a new sta- 
tionary state at any finite position x. A physical realization would be an abrupt change 
in the bias voltage of an electronic device. The characteristic times of the transient are of 
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practical interest to determine the transport properties of small mesoscopic structures, 
but modelling the process by means of a grid discretization of space in a "finite box" is far 
from simple ^230j , since the boundary conditions at the box edges are not known a priori 
and involve simultaneous injection to and absorption from the simulation (box) region. 
Some approximate ways to deal with the transients have been proposed ! 231|230|232|233j . 
An exact numerical solution generally requires non-local in time "transparent" bound- 
ary conditions [234j . Some cases are solvable analytically in terms of known functions: 
Delgado et al. [64] worked out an explicit and exact solution of the transition between 
stationary monochromatic waves due to an abrupt potential switch for a step potential 
~VoQ{x) that changes the step height suddenly to —VqQ{x), Vq > Vq > 0, and other po- 
tential profiles, e.g., containing square single or double barriers can be treated similarly. 
The basic trick to find the exact solution is to implement the action of the evolution 
operator of the new Hamiltonian on the initial state using an integral expression in the 
complex momentum plane obtained by Hammer, Weber and Zidell [235]. For the step 
barrier that changes its height abruptly one starts with a stationary wave function of 
the original Hamiltonian with incident (momentum qo), reflected and transmitted waves, 
ip{t = 0) = il'i+ip2 + '>p3- Each part is treated separately as a modified Moshinsky shutter 
problem and combined later. Let us analyze, for example, the incident case: The eigen- 
states of the final Hamiltonian can be written as the incident plane wave (for all x) plus 
a rest. The time evolution for ijji is then given by 



where q is the momentum from the left potential, Ci goes from — cx) to oo above all 
singularities and 4>i{q,t = 0) is the momentum representation of 'ijji{x,t = 0). From here 
the integral is solved by deforming the contour along the steepest descent path. This 
provides w-function terms and branch-cut contributions, but the latter can be neglected 
in many cases. The characteristic time scales for the transitions in the step potential admit 
simple classical analogs. Essentially the same method can be applied to the outer parts 
of the wave for more complex cut-off potentials such as a double barrier with a well that 
can be adjusted to make the wave on or off-resonance |65|201j . The contour deformation 
speeds up convergence and provides physical insight as discussed in 17.11 However the 
method is difficult to apply for the evolution of the wave component initially inside the 
quantum well region which was calculated with the systematic pole expansion of Garci'a- 
Calderon and coworkers [65j , see Section[6l In summary, the build-up process is dominated 
by the interference between the incident and resonance poles, leading to a non-exponential 
growth in which the resonance lifetime is a relevant characteristic time; the decay process 
is dominated by the displaced resonance; and the trapping process involves very long 
time transients due to interferences between a non-exponential resonance contribution 
and a bound- virtual state pair. Although a real quantum dot structure is much more 
complicated, a simple model potential, with realistic parameters, shows the importance 
of transient behavior over a wide (on the pico-second to nano-second and even longer) 
time scale. This work was extended to the case of interacting electrons at the mean field 
level by coupling Poisson and Schrodinger equations [66j . A critical coupling was found 
(determined by the electron sheet density in the well) above which permanent trapping 
does not occur. 




(132) 
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Julve and Um'es have similarly described the transient from a pure plane wave on 
all space and the corresponding state for a square barrier using Laplace transform tech- 
niques [236] , and Moshinsky and Sadurni have studied the dynamics of a deuteron sud- 
denly subjected to an electrostatic field [237j . Finally, Abdallah and Pinaud |238j have 
determined non-local in time "transparent boundary conditions" necessary to solve the 
coupled Schrodinger-Poisson equation for the transient evolution of a resonant tunnel- 
ing diode when the incoming amplitude representing the electron injection is prescribed, 
and Antoine et al. have recently reviewed different techniques to solve numerically the 
time-dependent linear and non- linear Schrodinger equation on unbounded domains [234j . 

9. Other transient phenomena 

9.1. Ultrafast propagation and simultaneous arrival of information in absorbing media 

Delgado, Ruschhaupt and Muga have shown that the temporal peak of a quantum 
wave may arrive at different locations simultaneously in an absorbing medium 1361221) . 
The arrival occurs at the lifetime of the particle in the medium from the instant when 
a point source with a sharp onset is turned on, or when a Moshinsky shutter is opened. 
(The effect remains for smoothed pulses [221 .) Other characteristic times may be iden- 
tified. In particular, the "traversal" or "Biittiker-Landauer" time (which grows linearly 
with the distance to the source) for the Hermitian, non-absorbing case is substituted by 
several characteristic quantities in the absorbing case. The simultaneous arrival due to 
absorption, unlike the Hartman effect, occurs for carrier frequencies under or above the 
cut-off, and for arbitrarily large distances. It holds also in a relativistic generalization 
but limited by causality, i.e., limited to distances where the time of arrival of the peak 
is larger than the time of the very first front. The ubiquitous peak is dominated by the 
saddle-point contributions above the cutoff frequency. A proposal for an experimental 
verification is described in |36| making use of a quantum optical measurement of atomic 
densities in a metastable state by fluorescence. Another interesting possibility is the use 
of absorbing wave guides and electromagnetic waves. One advantage of the waveguide 
realization with respect to quantum particles is that the effect could be measured in a 
non-invasive way. It also makes it possible to trigger several devices at the same time or 
transmit information so that it arrives simultaneously at unknown locations. This cannot 
be achieved by standard transmission methods because a receiver could not resend an 
information bit to the closest one faster than the velocity of light in vacuum c; nor can we 
design the timing of a series of signals from the source so that they arrive simultaneously 
at different receivers if their locations are unknown. 

The simplest case is that of a source with a sharp onset, and the treatment is very 
similar to the one used in l7.2i substituting the real potential by a complex one, 

fcg becomes also a complex number so that the role of tbl in the formal solutions (see 
Eq. (|120p ) is now played by the complex number 

TcBL = mx/{-ikoh), (134) 
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and the steepest descent path in the fc-plane crosses the pole in fco at a time Tc, Eq. 
(|114p . For carrier frequencies below the cut-ofF the saddle contribution dominates up 
to exponentially large times as in the Hermitian case [34j . so that Tc is not of much 
significance whereas, above threshold, gives approximately the arrival of the main 
peak. From d\%lj(x, t)\'^ jdx = we may obtain the position x{t) of the "spatial" maximum 
at time t. This also defines (by inverting x{t)) a function ts{x)^ namely, the time when 
this spatial maximum arrives at x. One finds that the spatial maximum in the large- a; 
region is given by ts — \tcbl\i a role played by the real tbl in the absence of absorption 
[239) . 

Thus a single real quantity tbl hi the Hermitian case has been substituted, for a 
non-zero V/, by three different quantities: the time for pole-cutting Tc, a complex Tcbl, 
and its modulus |tcbl|, all of which tend to the "Biittiker-Landauer" traversal time 
Tbl = TcBhiYi — 0) without absorption 189^ 

Similarly, we may fix x and solve d\'ip{x^t)\'^ / dt = to get Tt{x), namely, the time of 
the first "temporal" maximum. At variance with the Hermitian case, the time of arrival 
of the maximum is not proportional to tbl {ti = tblI^/^ when Vi =0). The equation 
9|?Asp/9i = cannot be solved analytically for V/ ^ in a generic case, so there is no 
explicit formula for Tt{x). Nevertheless, if x » 2~^/^|fco|/i^/mVf , one finds 

n « A, (135) 

i.e., the temporal maximum coincides with the mean "life time" hj^Vx of a particle 
immersed in the absorbing potential, and it is independent of x and wq, which is the most 
important result of [36j . A time- frequency analysis [217134] shows that the maximum is 
not tunneling but it is dominated by frequencies above the cut-ofF, in particular by the 
saddle frequency corresponding to the "classical" velocity required to arrive at x at time 

Tt(x). 



9.2. Breakdown of classical conservation of energy in quantum collisions 

In most quantum interferometers the phase differential between paths that meet in a 
coordinate-space point or region at a given time lead to constructive or destructive wave 
combinations and thus to fringes, but the paths can also interfere in momentum space 
and produce momentum-space fringes. An example of interference in momentum repre- 
sentation is the collisional, transitory enhancement of the high momentum components 
of a wave packet |240|241|199|202j . This enhancement leads to a violation of the classical 
conservation of energy since the probability of finding the particle with a momentum 
larger than a given value may exceed the classical bound during the collision. The effect 
was shown to be quantitatively significant for incident energies well above the barrier top 
when a transient suppression of the main incident momentum components "pushes" the 
momentum distribution towards lower and higher values 199 202 . This was explained 
as an interference in momentum representation between components associated with in- 
cidence and transmission. A simple picture emerges by assuming that in the midst of 
the collision the wavepacket density is unaltered with respect to free motion, but the 
transmitted part acquires the phase of the transmission amplitude e*"^^ . Adjusting this 
phase with the incident or barrier energies constructive or destructive interferences ap- 
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pear in momentum space. For Gaussian states impinging on square barriers, an analytical 
approximate expression for the wave function is [199] 



ip{k, t) oc w[ui{k, t)] + T{k)w[—UT{k, t)] 



(136) 



where uj{k,t) and UT{k,t) are known functions of momentum and time [199] . and T(fc) 
is the transmission amplitude. If, for the central momentum components, the modulus 
|r| is close to unity and its phase is an odd multiple of tt, the two contributions in Eq. 
(jl36p interfere destructively to cause the transitory suppression of the main momentum 
components of the incident wave packet. This is rather surprising from a classical per- 
spective since only a negligible part of the wave packet is located at the potential barrier 
or, equivalently, most of the wave packet is on regions with zero potential energy. Both 
transmitted and incident components are dominated by the central momentum of the 
distribution, but instead of adding up classically, they interfere destructively. 

By sending the wave packet towards the appropriate barrier potential and then switch- 
ing the potential off at the time when the maximum interference occurs, this effect may 
be used to convert a Gaussian-shaped wave packet into a bimodal one in momentum 
representation, since the free evolution afterwards will not modify the momentum dis- 
tribution. However, the transitory character of the effect could hinder its observation or 
its use as a preparation tool for controlling the atomic state of motion in atom optics 
or quantum information applications. To remedy this problem an extra potential barrier 
was added in 202 . The interference effect becomes more complex but it is possible to 
increase its duration significantly making it more robust and easier to observe. 

For an alternative to the scattering configuration using trapped atoms and "phase im- 
printing" by detuned lasers which illuminate only half the atomic cloud see |242j . If a tt 
phase is imprinted, the interference dip formed in momentum space is centered at fc = 0, 
but the minimum is displaced linearly if the phase deviates slightly from tt. Turning off 
the external trap thus creates a dark notch in coordinate space moving with controllable 
velocity, similar to the "dark solitons" observed experimentally ^243 . In j242j the pertur- 
bation of imprinting imperfections and many body effects are investigated as a necessary 
step to determine the potential applicability of momentum space interferometry. The 
imprinted phase carries information about the laser interaction (time, laser intensity, 
frequency) that can be obtained from the observable notch. 

9.3. Effect of classically forbidden momenta 

Suppose that a classical ensemble of independent particles in one dimension is initially 
confined (at t — 0) in the spatial interval a < a; < 6 < 0, and allowed to move freely after 
t = 0. Only particles with positive momenta may arrive at positive positions for i > 0. In 
contrast, the quantum wave function involves negative- momentum contributions as well. 



oo 




(137) 



— oo 



where Ek ~ k^ti^ /2m, and 
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V'(fc) (27r)-i/2 / dxe-'''''^P{x,0) (138) 



is the wavenumber representation of the initial state. The effect of negative momenta for 
a;, < > is however a transient one; the total final probability to find the particle at a; > 
is given only by positive momentum components, 

Pr(oo)=hm j dx\i^{x,t)\^ ^ Jdk\ip{k)\^, (139) 



see e.g. [244) or [245) . since there are no bound states. This negative-momentum effect 
is also present in collisions, where it may be combined with other classically forbidden 
effects [246] . Consider the family of cut-off potentials of the form 

r if a; < c , 

V{x)^lu{x), if c<x<d, (140) 
I Vb, if a; > d, 

where Vq > and V{x) are real. Let us suppose that the maximum value of the potential 
is Vm- For a classical ensemble confined between a and b (such that a<b<c<d<0) 
only particles with initial momentum above the barrier, k > kM, km = (2mVM)^^^/^, 
may pass to the right of the potential region. In the quantum case however, there are con- 
tributions from all the components of the initial wave packet: (a) k > kM = i'2mVMy^^/fi 
(above-the-barrier transmission); (b) kg < k < k^.j ("asymptotic" tunneling, these mo- 
menta contribute to the transmission probability Pt (oo)); (c) < k < ko ("transient 
tunneling", these momenta do not contribute to Pt(oo)); (d) —oo < k < (transient 
negative-momentum effect); (e) kj — ijj, 7j > (tunnel effect associated with bound 
states, which do contribute to Pt{oo)). 

The negative momentum effect is the only one that survives for free motion, and it has 
been frequently overlooked, exceptions being [235l247l248l246lll9j . 

9.4. Time of arrival and Zeno effect 

An elementary model for the time of arrival measurement at a point in space, a; = 0, 
consists in performing instantaneous and frequent projections [245j ipi^: = i^i^j tj -)^{~ 
with the rate of norm lost playing the role of an operational time-of-arrival distribution. 
After each chopping, the time evolution can be expressed as a combination of Moshinsky 
functions Jj)7j, and from their short-time asymptotics it follows that when the projec- 
tions are performed very frequently the wave packet is totally reflected, according to the 
quantum Zeno effect. The first discussions of the Zeno effect, understood as the hindered 
passage of the system between orthogonal subspaces because of frequent instantaneous 
measurements, emphasized its problematic status and regarded it as a failure to simulate 
or define quantum passage-time distributions |245l249j . Muga et al. however [250], have 
recently showed that in fact there is a "bright side" of the effect: by normalizing the little 
bits of norm N removed at each projection step, 

Ilzeno = Hm — ^^(^ , (141) 

st^o 1 - N{oo) 
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a physical time distribution defined for the freely moving system emerges, which turns 
out to be in correspondence with a kinetic energy density j251| . rather than with a 
time-of- arrival distribution. 

10. Experiments 

Compared to the abundant theoretical work on elementary quantum transients, the 
pace of experimental observations has been relatively slow, but fortunately the trend is 
changing very rapidly. One reason for the difficulties is the feebleness of diffraction in time 
in some systems. Dissipation [38], environmental noise |15j . and repulsive interatomic 
interactions [IB] tend to suppress DIT-like phenomena. Moreover, the smoothing of the 
aperture function of the shutter [12| leads to apodization in time, which delays the 
manifestation of quantum transients until the revival time [15' discussed in Section [4] 
Similarly, DIT may be suppressed by superposition of carrier frequencies. A finite width 
Av of the velocity distribution induces a spreading Avt in coordinate space which should 
be negligible compared to the width of the main fringes in ([5]). Hence, DIT weakens 
with the width- broadening of the velocity distribution [48l34j . which often arises as a 
direct consequence of the uncertainty principle. In particular, it comes into play in the 
dynamics of beams of finite-size and velocity selection processes |7I8] as well as localized 
states under confinement ^9)10|48j . 

The analysis of Goldemberg and Nussenzveig [7] rightly pointed out that the experi- 
mental observation of quantum transients would be simpler for the time-energy uncer- 
tainty principle than for interference effects on the density profile, though the latter can 
be enhanced in matter-waves with attractive interactions |152j and using constructive 
interference with reflected components from moving mirrors |42j . 

Despite the hindrance, a number of remarkable experiments have focused on quantum 
transients and time-dependent optics of different matter waves. A series of experiments 
[26127128] culminated in the confirmation of the time-energy uncertainty relation derived 
by Moshinsky for pulse formation 24J. The matter wave source consisted of caesium 
atoms in a magneto-optical trap (MOT) which are released and fall under gravity to 
bounce off of a mirror, made of an evanescent laser wave. A selection of atoms with 
a given total energy was achieved in a first bounce, while the double-slit in time was 
implemented during the second bounce, after which the resulting pulse was probed. In 
such a way the energy distribution of the beam when chopped in pulses was measured in 
agreement with the time-energy uncertainty relation 23 24J. Moreover, it was possible to 
register the interference pattern in the arrival-time distribution of the cold atoms after 
reflection at the double temporal slit. A similar setup was used by Colombe et al. [75] in 
which a Bose-Einstein condensate, with negligible non-linear interactions, was dropped 
over a harmonically vibrating mirror which can act as a phase modulator [26] . After the 
bounce the atomic carrier and sidebands were directly measured by absorption imaging 
of the density proflle [75] . 

Cold atoms imply long wavelengths and time scales favorable to the observation of 
transients. The recent development of guided atom lasers [2531254] , may pave the way 
for the direct observation of the DIT fringes of a suddenly turned on matter wave beam 
[151121] , allowing for a close implementation of the initial setup discussed by Moshinsky 

i- 
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Historically, ultracold neutrons were initially the favorite matter waves to study such 
phenomena. In this case, the interaction with a potential periodically modulated in time, 
was shown to induce frequency side-bands. The resulting energy spread of the matter 
wave could be detected using a high-resolution time-of flight instrument [71|21I72] and 
was shown to be in agreement with previous theoretical insight [68|69|70|13ll4j . The 
interaction with a grating moving across a beam of ultracold neutrons similarly induces 
an energy spreading measured in [252] . 

Further studies on the time-energy uncertainly relation have been performed at the 
attosecond scale during the photoionization of Argon atoms. A double slit in time was 
implemented using few-cycle ionizing fields, where the momentum distribution of the 
resulting photoelectron carries the information of the time at which the ionization occurs. 
The interference in the energy distribution was then measured as a function of the time 
interval between the generation of different isoenergetic photoelectrons |73| . 

A boost in the experimental progress on quantum transients has recently taken place 
[76j exploiting the analogy between paraxial optics and the time-dependent Schrodinger 
equation |255j . By means of this analogy, the bifurcation in the density profile of freely 
expanding excited states of a hard-wall trap predicted in [15] (see section [3]) has been 
reported. The free propagation of eigenstates and coherent states of square quantum bil- 
liards (essentially a 2D hard- wall potential), was simulated using the transverse modes of 
oxide-confined vertical-cavity surface-emitting lasers and recorded with a CCD camera. 
The analogy with matter- wave dynamics holds under the replacements t — > z, where the 
z-axis is the direction of vertical emission, and m/h 27r/A, A being the lasing wave- 
length. Moreover, for a laser cavity chaotically shaped a random branching interference 
of the coherent lasing modes was observed. 

One more recent new route for experimental access to quantum transients and DIT 
phenomena is based on the time evolution of the population of an atomic level, i.e., 
transients do not only occur for translational degrees of freedom but also for internal 
degrees of freedom after or during laser excitations. Specifically, linearly chirped pulses 
where the instantaneous frequency changes in time passing trough the resonance con- 
dition for a particular transition of rubidium atoms have been shown to produce DIT 
oscillations ( "coherent transients" ) of the excited state in the low field regime [77|256j . 
They are probed with pulses by a second laser on the subpicosecond time scale and can 
be attributed to the interference between resonant and off-resonant excitations. 

These later developments and, in addition, the fact that techniques of ultra-short laser 
pulses in the attosecond time scale (without the need for high-intensity X-rays) are a 
gateway to a plethora of time-dependent electronic processes that were too fast to be 
observed, [257] make the prospects for quantum transient experiments rather brilliant. 

11. Final comments 

More than fifty years after Moshinsky's pioneering investigation on the quantum shut- 
ter, quantum transients keep posing theoretical and experimental challenges and fasci- 
nating open questions. The initial interest in nuclear physics has shifted to applications 
in cold atoms, quantum optics, and semiconductor structures. Recent experimental capa- 
bilities to prepare and control light pulses, atoms, electrons and their interactions make 
possible a whole new world of quantum transients that has been barely explored and 
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will provide insight into many aspects of fundamental physics that had been until now- 
difficult to probe or change. Few and many-body physics will offer, for example, ample op- 
portunities for discovering and exploiting interesting phenomena without single-particle 
counterparts. At the single particle level, new techniques, ideas, and physical realiza- 
tions of quantum transients set a sound ground for realizing many effects described in 
the theoretical studies, and will itself motivate further theoretical work. A challenge is 
the direct observation of matter-wave diffraction-in-time oscillations, but promising and 
versatile substitutes are the electromagnetic analogy that has been experimentally real- 
ized near the completion of this work, and "coherent transients" of atomic populations 
under controlled field excitation. Applications such as the control and understanding of 
quantum dot dynamics, quantum gates operation, tunneling, atom lasc^rs. momentum in- 
terfcromctry, coherent control, or simultaneous transmission of information hold promise 
of a vigorous development of research on quantum transients, technological advances and 
exploration of new physics. 
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Appendix A. Propagators in quantum transients 

Many relevant quantum transients can be formulated as an initial value problem. 
A rather general setting considers a given eigenstate \t/j{t = 0)) of a Hamiltonian Hq 
which is suddenly perturbed to a new Hamiltonian H. The initial state undergoes a 
non-trivial evolution \il>{t)), which for non-relativistic particles is governed by the time- 
dependent Schrodinger equation. Formally, the time-dependent state vector becomes 
\tp{t)) = U{t, 0)\t{j{0)) where the evolution operator for Hermitian Hamiltonians satis- 
fies unitarity, WU = UW = 1, and the composition property U{t,t')U{t' ,0) = U{t,0). 

The most general form of the time-evolution operator is given by the Dyson series, 

U{t, t') = l + Y^ (^1^ Jdtijdt2--- J dtnH{ti)H{t2) . . . H{tn), (A.l) 

t' t' t' 

though as it is generally the case specific problems admit simpler forms. Indeed, for a 
time-independent Hamiltonian, U = exTp{—iHt/h); while, iidtH but the Hamiltonian 
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at different times commute {[H{t), H{t')] = for any t and t'), U(t, 0) = exp(-i dt'H(t')/h). 
In coordinate representation, and using the resolution of the identity, 

{x\m) - / dx'{x\U{t,t')\x'){x'm')), (A.2) 



one naturally arrives at the concept of the propagator, which is the kernel K{x^ t\x' , t') = 
{x\U{t, t')\x'), this is, the coordinate representation of U. In the general case in which the 
Hamiltonian has both a discrete set of eigenstates labelled by the quantum number 
n, and scattering states {ve} of energy E, the spectral decomposition of the propagator 
reads 

K{x, t\x', t') = J2 04a;)0„(a;')*e-^^"(*-*')/'^ 

n 

oo 

+ / dE^Eix)ipE{x're''''^'-'">/\ (A.3) 



where * denotes complex conjugation. A path integral representation of the propagator 
is possible, 

A-,Ml<.0 = /l'.[.lexp(i5H), ,A.4, 

in terms of the classical action S[x] — J^, dT{mx^/2 — V{x)). The propagator is thus an 
integral over all paths x[t] satisfying the boundary conditions x[t'] = x' and x[t] ~ x. 
Without intention of being exhaustive, we next point out some specific features relevant 
to the study of quantum transients: 

- The propagator for quadratic Hamiltonians acquires the simple form '258] 

K{x,t\x',t') ^A{t)exp(^^Sc.i[x]^ , (A.5) 

where 5*2; [x] is the action corresponding to the classical path, as it follows from the 
stationary phase approximation. Here A(t) is just a time-dependent function indepen- 
dent of the position variables. For bounded regions, other propagators are known to 
collapse into a countable sum of classical paths [259] . 

- Hamiltonians with an arbitrary time-dependent linear potential can be mapped to the 
free particle problem [260] . 

- Perturbations of the delta function type can be taken into account via path integral 
summation of perturbation series ^261j and the Laplace method (110117] . 

- There are important classes of time-dependent problems which can be mapped to time- 
independent ones making use of integrals of motion or the so-called Duru's method 

- The propagator for a problem with Dirichlet, von Neumann and mixed boundary 
conditions can be obtained from the unconstrained one using the method of images 
[258 263 264 261 265 . Simple examples involve a hard- wall at the origin, K^aii [x, t\x' , t') 
K{x, t\x' ,t') — K{—x, t\x' , t'), where K is for instance the propagator for the free Hamil- 
tonian or harmonic oscillator. Similarly, the dynamics of a particle in a box can be 



59 



described by taking into account infinite images, even in the presence of moving bound- 
aries |266I261| . 

The general theory for the path integral approach to quantum mechanics can be found in 
[258|263|264j . In particular, explicit expressions for propagators are collected in [261j . an 
excellent account of exactly solvable path integral problems. A recent extension of some 
of these results exploits methods of supersymmetric quantum mechanics |267|268|269) . 
Finally, one should mention representations of the retarded time-dependent Green func- 
tion as expansions involving resonant states [2801287] , which are summarized in Appendix 

m 



Appendix B. The Moshinsky function 



The Moshinsky function arises in most of the problems where ID quantum dynamics 
involves sharp boundaries well in time or space domains |4I5I24I16] . Indeed, it can be 
considered "the basic propagator of the Schrodinger transient mode" |274j . Similarly, it is 
found when considering free propagators perturbed with point-interactions. Its standard 
definition reads 



M{x, k,t) 



e 2ht 



w{-z), 



(B.l) 



where 



1 + i ht 



ht ) 

An equivalent expression is 



M{x,k,t) 



-erfc 



X — hkt/m 
m 



(B.2) 



(B.3) 



Its absolute square value, \M{x, k, t)p, admits a simple geometric interpretation in terms 
of the Cornu spiral or clothoid, which is the curve that results from a parametric represen- 
tation of the Fresnel integrals. Indeed, another frequent representation of the Moshinsky 
makes use of the complex Fresnel function T , 



M{x,k,t) = -j=e-'^+'''^ 



(B.4) 



where J-{9) = exp^inu^ /2)du — C{9) + iS{9) is an odd function of its argument 



k 



mx 



mn \ ht J 

Directly from Eq. (|B.4p . or using Eqs. (|C.6p . (|C.7p . and (jB.ip . it follows that 



|M(x,fc,i)|' 



1 



1 2 



S{9) 



(B.5) 



(B.6) 
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which, up to the definition of 9, is identical to the resuh from Fresnel diffraction from a 
straight edge [19120] . It is possible to understand physically the Moshinsky function as a 
freely evolved cut-off plane wave. Using the free propagator Q, 







Mix,k,t)^ I dx'Koix,t\x',0)e''''''ei-x')^ I rfxJ-^e' "'^^^*''' +'^"'. (B.7) 
An integral representation in the fc-variable, follows from the same physical problem, 

- p ^—ihk t/2m-^ik x j C p~'^^ — -+ik x 

-OO r+ 

where Tj^ is a contour in the complex z-plane which goes from — oo to cx3 passing above 
the pole. 

Many properties of the Moshinsky function directly follow from its relation to the 
Faddeyeva w{z) function. In particular, under inversion of both x and fc, it satisfies 

M{x, k, t) + M{-x, -k, t) = 6*'=^-*^*. (B.9) 

Similarly, the asymptotic behavior of the Moshinsky functions for \x — hkt/m\ — > oo 
can be found from those of w(z) when z oo |34j . For x < hkt/m, and z as in Eq. (|B.2p . 
one finds in terms of the Gamma function T{y) that 



Mix,k,t)^e'^'=~'— + 



2iy/TTZ 



1 + E 



(2n- 1)!! 



^ikx—i 



^ (2z^) 

n=l ^ ' 

Hk^t e'^R*" 1. - r (n 



2\n 



2711 

ri=0 



whereas in the complementary region x > hkt/m only the series survives. 

B.l. Integral and differential equations 

The following relations are known [110] 

d^M{x,k,t)^ikM{x,k,t) ~ Ko{x,t\0,0), (B.ll) 

dkM(x,k,t) ^ikM(x,k,t) - —dxM(x,k,t), (B.12) 

m 

dtM{x,k,t)^^dlM{x,k,t), (B.13) 
2m 

= Y^Koix, t\0, 0) ~ ^d^Mix, fc, t), (B.14) 



where Ko{x,t\0^0) — yiSl^* ' ^^^'^ propagator from the origin. In ad- 
dition, the indefinite integral (which arises in tunneling problem through short-range 
potentials |110|17j ) holds 
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M{ax' + 6, c, t)e"'^ da; = t^^^ y [M(aa; + 6, c, i) - M(aa; + h, -q, t)] , 

(B.15) 

but for q = —ac, when I'Hopital's rule is to be applied. 
Appendix C. The to-function 

The Faddeyeva function is also known as the complex error or plasma function [78] . It 
is an an entire function related to the complementary error function [79], 

w(z) e"^'erfc(-iz). (C.l) 

Moreover, it admits the integral representation 

1 f e""^ 

w{z) = — / du , (C.2) 
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where r_ is a contour in the complex z-plane which goes from —oo to oo passing below 
the pole. Actually, it satisfies the following equation, 

oo ^ 

/ du = i7rsign(3z)ii;[sign(3z)2:], (C3) 

J u z 

— oo 

and it obeys the so-called Faddeyeva identity, which follows from Cauchy's theorem 

+ w(-z) = 2e"^'. (C.4) 

Moreover, 

w{z*) ^w{~z)*. (C.5) 
We further notice that the Frcsncl integrals 



S{9) = / cos —dt, 



C{9) = I cos— dt, (C.6) 



are related to the w-function by [75], 



It has a series expansion 



1 _ e-^V2^ ( i±l^i/2^ 



(C.7) 



n=0 ^ \2 ^ ^) 
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whereas, its asymptotic expansion reads (expanding {u — z) ^ around the origin in Eq. 
(|C.2p and integrating term by term) 



w{z) 



1 + E 



{2m- 1)!! 



- (2z 



2\m. 



2e"^ e(-3z), (C.9) 



for z oo. We also note that its derivatives can be evaluated in a recursive way, using 
the relation 

w'{z)^~2zw{z) + ^. (C.IO) 

In the upper half-plane, < |w(2:)| < 1, and its efficient computation is described in 
[270|271j . The real part of the w-function is proportional to the Voigt spectral- line profile, 
a convolution of a Lorentz and a Gaussian profile. 

Apart from its close connection with the Moshinsky function, w-functions arise also 
naturally in the dynamics of some initial states which are not truncated in position 
space: for Gaussian- like examples in free motion see |272| . and for Lorentzian states in 
momentum representation [llOj impinging on a separable interaction see [273) . 

Appendix D. Inhomogeneous quantum sources 

The time-dependent Schrodinger equation for a freely moving particle with an inho- 
mogeneous source term switched on at i = reads 

^ = ^^a^+2;;^a^- 
Using the free Green's function, 

Go(x, t\x\ t') = Q{t - t')Ko{x, t\x't'), (D.2) 

which obeys 

LGoix, t\x',t') = ihS{x - x')6{t - t'), (D.3) 
the solution for 'ip{x, t > 0), assuming 'ip{x, i < 0) = 0, is given by 
t 

ih 



il;{x,t) = ^ I dt' I dx'Go{x,t\x',t')a{x',t'), (D.4) 



as can be verified by substitution, note the upper limit t+ in the time integral. 

In what follows we shall be interested in the dynamics of an inhomogeneous source of 
the form a{x, t) = e~'"*£i(x), where hui = (Sfc)^/2m. This can be related to the Moshinsky 
shutter problem with the initial condition ip{x,0) = sm(kx)Q{—x), which evolves into 
ip{x, t) = [M{x, k, t) - M{x, -fc, t)]/{2i). We note that 

^{x,t) = ^ f g{q)e^'^^-^'^dq, (D.5) 



/2tt 



63 



where 



5(g) = 



/2^[(q + ie)2 - fc2] • 



(D.6) 



has two poles in the lower half-plane, close to the real axis. The integral can be deformed 
along the contour r+ which goes from —oo to oo passing above the poles at q = ±fc, 



e''"'dq. 



(D.7) 



Rewriting 



e 2" 
q^ — fc2 



2, 



dt' 



e 2m 

^2 _ ^2 ■ 



only the first term contributes along the contour r+, and therefore 



hk 

2im 



1 

2^ 



- hq^ it-t' ) 



(D.9) 



Identifying the term in brackets with the free Green's function and comparing (|D.9|1 with 
(jD.4p . it follows that 



h^k 



(D.IO) 



In other words, the DIT setup with a totally reflecting shutter (reflection amplitude 
R = —1) is equivalent to an inhomogeneous point-like source suddenly turned on. This 
result was first obtained by Moshinsky |216j . 



Appendix E. Resonant states 



The first ideas for a theory of resonant states originated at the end of the twenties 
of the last century with the work of Gamow on nuclear radioactive decay |275j . He 
considered a situation where a particle, initially confined inside a three dimensional 
region by a potential barrier, goes out of it by tunneling. In order to describe the above 
process, Gamow restricted the discussion to spherically symmetric systems and looked for 
solutions to the Schrodinger equation which at large distances consist of purely outgoing 
waves. He realized that the absence of incoming waves in the solution is fulfilled if the 
energy eigenvalues are complex. In 1939, Siegert used these purely outgoing boundary 
conditions to derive a dispersion formula for elastic scattering by a potential of finite range 
|276j . The same definition of resonant states was also used by Humblet and Rosenfeld 
to formulate a theory of nuclear reactions [277j . One sees that resonant states provide a 
unified description of decay and scattering problems. Following an idea by Peierls [278], 
subsequent developments involved a consideration of the analytical properties of the 
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outgoing Green function to the problem |279I280| . These led to expansions of this function 
and to a consistent normalization condition for resonant states |280j . The extension of 
resonant states to one dimension (i.e. to the full line) was made in Ref. |281| . Along 
the full line, the fact that there are two edges to the potential leads to a modification 
of expression for the normalization condition |282|61j . For the half- line the resonance 
formalism is similar to the case of zero angular momentum in three dimensions. 

Here we present some properties of resonant states relevant for the discussion of sections 
El and [71 

Most of the work on resonant states in one dimension refers to interactions of finite 
range, i.e., extending along a finite region of space. A given resonant state satisfies the 
Schrodinger equation of the problem with complex energy eigenvalues [61j . 

-^Un{x)+[kl-U{x)]u„ix)=0 (E.l) 

where fc^ = 2mEn/h'^, with the complex energy En — — ir„/2, U{x) — [2mV{x)]^^'^ /H 
where V{x) vanishes beyond the interval < x < L. The solutions to the above equation 
satisfy outgoing boundary conditions at a: = and x = L, given respectively by 



and 



d 

dx 



d 

dx 



= -iknUniO), (E.2) 

x=0 



= iknUn{L). (E.3) 



The amplitude of a resonant state increases exponentially with distance beyond the in- 
teraction region, i.e., since fc„ = a„ — i/9„, say for x > L, Un{x) ~ exp(iQ;a;) exp(/3a;), 
and therefore the usual rules concerning normalization, orthogonality and completeness 
do not apply for resonant states. For that reason it was generally believed that these 
functions were not very useful for calculations. As mentioned above, developments in- 
volving the analytical properties of the outgoing Green function to the problem have led 
to a consistent normalization condition and completeness condition for resonant states. 
Notice also that at distances beyond the interaction region, say x > L, the resonant state 
may be written as 

Unix, t) - e^(o'^^-£^t/k)^P^x^-T^t/h^ (2 4) 

which represents an outgoing wave of wave number an associated with a state of positive 
energy En- The second factor shows that this state decays at the rate Tn/ti, in agreement 
with the usual interpretation of the width of a resonance level. Since the velocity of the 
outgoing particle is w„ — han/m, one sees that the exponential increase with x of the last 
factor in Eq. I|E.4[) . may be interpreted as due to the fact that at a distance x one finds 
those particles which left the interaction region at a time t — x/vn when the amplitude 
there was larger by the factor exp[r„a;/(2fiu„)] = exp(/3„a;). 

It is well known that for finite range interactions the outgoing Green function G^{r, r';k), 
as a function of k, can be extended analytically to the whole complex fc-plane where it 
possesses an infinite number of poles distributed in a well known manner [2831284] . Purely 
imaginary poles, situated on the upper half fc-plane, correspond to bound states of the 
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problem whereas purely imaginary poles seated on the lower half /c-plane are related to 
antibound states, also called virtual states by some authors. On the other hand, complex 
poles are only found on the lower half fc-plane and are associated with resonant states. To 
each complex pole at k„ = a„ — with a„, /?„ > there corresponds, from time rever- 
sal invariance, a complex pole fc_„ situated symmetrically with respect to the imaginary 
axis, i.e., fc„„ = Also u_„(a;) = "«* (a;). 



E.l. Time-independent resonance expansions 



The outgoing Green function to the problem satisfies the equation 



dx"^ 



G+{x,x';k) + [fc2 - U{x)\ G+{x,x';k) 



2m 



S{x — x'), 



(E.5) 



where fc^ = 2mE/h'^, with E the energy. The solution to Eq. (jE.Sp satisfies outgoing 
boundary conditions at .x = and x — L given respectively by 



and 



^G+{x,x';k) 
ox 



—G+{x,x';k) 



= ~ikG+{0,x';k) 



= ikG+iL,x';k). 



(E.6) 



(E.7) 



It may be shown that the residue Cn{x,x') of G^{x,x';k) at a pole fc„ seated on the 
complex fc-plane is given by [61j 



, / 2m\ Unix)Un{x') 

C^{x,x)-^j^j — , 
provided the resonant states are normalized according to the condition 

L 



ulXx)dx + i 



ul{0) + ul{L) 



1. 



(E.9) 



Notice that for a bound state (i.e., fc„ = i7„) the contribution of the two terms on right- 
hand side of Eq. (|E.8|) corresponds exactly to two integral terms that allow to express 
the normalization condition for the bound state in the usual form u^{x)dx = 1. 

An interesting expression follows by considering the Green theorem between Eq. (jE.ip 
for Un{x) and a similar equation for Um(x) and then using the corresponding boundary 
conditions for Un(x), i.e., Eqs. (jE.2p and (|E.3p . and similarly for Um{x), to obtain 



, ^ , , .Un{0)Um{0) + Un{L)Um{L) 
Un[X)Um[X)dx + I , , = 0. 



(E.IO) 



The above results suggest to look for an expansion of the outgoing Green function of 
the problem in terms of the set of resonant states. This may be achieved by considering 
the following integral. 
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1 



G+ix,x';k) 



dk' , 



(E.ll) 



where C represents a large closed contour of radius L in the fc-plane about the origin in 
the clockwise direction which excludes all the poles, fc„ and the value at k' = k. Using 
Cauchy theorem it follows that I — and hence one may write 



1 



G+{k') 
k' -k 



dk' 



Cr 



= 0, 



(E.12) 



where for simplicity of notation G+(x,x'; k) is denoted by G^{k') and Cr represents a 
large circle centered at the origin; the contours C„ encircle each of the poles fc„ that are 
enclosed by Cr, and the contour Ck encloses the value k' = k. All these contours are in 
the counterclockwise direction. 

Using the theorem of residues and Eq. (jE.Sp one may write Eq. (IE.12P as 



G+{x,x']k) = 



2m 



N 

E 



Un{x)Un{x') ^ J_ 



G+{x,x';k) 
k' - k 



dk'. 



(E.13) 



Ca 



The above equation is defined for all x,x' > 0. However, for many applications we need 
to know G'^{x,x';k) only along the internal region of the interaction. In this case it 
has been proved rigorously along the half-line that the outgoing Green function goes 
exponentially to zero along all directions in the complex /c-planc for all values of x and 
x' along the internal interaction region, except at the points x = x' = ~ L |285l286j . 
This is denoted as ^ {x, x') ^ L i.e., 



G'^{x,x';k) as |fc|— >oo. 



^ {x, x') ^ L. 



(E.14) 



Therefore by extending the radius R of the contour Cr up to infinity one then obtains 
that the integral term along Cr^ in Eq. (|E.12p vanishes exactly. Hence, one is left with 
an infinite sum over the full set of resonant terms, 



2m\ 



u„{x)u{x'„) 



^ (x, x') ^ L. 



n— — oo 



2kji(^k kn) 

Substitution of Eq. (jE.lSP into Eq. (IE.5P implies the fulfillment of the relations 
g u^ixWix') ^ ^ ^ 

n— — oo 

and 

1 °° 

- ^ Un{x)Un(x') ^ 5{X - x'), Q<{x,x')<L. 



(E.15) 



(E.16) 



(E.17) 



The first of the above two expressions indicates that resonance states are not independent 
of each other and the second one may be seen as a modified closure relationship. 
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Let us now consider a relationship between the outgoing Green function and the con- 
tinuum wave function of the problem. The continuum wave function satisfies the equation 

(P 

-TT^ik, x) + [fc2 - [/(x)]V'(fc, x) = 0. (E.18) 

For a particle approaching the potential from the left {x < 0) , the solutions to the above 
equation read, respectively for x < and x > L, 

= e^*="+i?(fc)e-^'=", (E.19) 

and 

=r(fc)e*'=^, (E.20) 

where R{k) and T{k) are, respectively, the reflection and transmission coefficients to the 
problem. It then follows by applying Green theorem between Eqs. (jE.lSP and (jE.Sp . and 
using the boundary relations given by Eqs. (|E.2p and (jE.3[) . that 

M^x) = ( —\2ikG+{Q,x]k), 0<x<L. (E.21) 

\ 2m / 

By considering the solution of the wave function at x = L, given by Eq. (|E.20p . into the 
above expression yields for the transmission amplitude 

T(k) = ( —] 2ikG+(0,L;k)e-'''^, (E.22) 
\2m / 

and similarly, using Eq. (jE.19[) . yields for the reflection amplitude 

R{k)^ 2ikG+ (0,0; k)-l. (E.23) 

\2m J 

Using Eq. (|E.15p into Eqs. (|E.2ip and (|E.22[) then leads, respectively, to resonance ex- 
pansions for the wave solution and the transmission amplitude, namely, 

oo 

%l){k,x)^ik ^ CnUn{x) 0<x<L, (E.24) 

n— — OO 

where Cn = Un{0)/[kn{k - kn)], and 

nk)-^k ± ^ff^e--. (E.25) 

n— — oo ^ ' 

The expansion for the reflection amphtude requires of subtraction terms. For the trans- 
mission amphtude one may alternatively expand G'^{x^x';k)exp{—ikL)^ to obtain 

nk)=^k ± ^ff^e--"^ (E.26) 

n— — oo ^ ' 
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E.2. Resonance expansion of the retarded time- dependent Green function 

The solution \^{t)) to the time-dependent Schrodinger equation may be written in 
terms of the retarded time evolution operator of the problem exp{—iHt), where t > 0, 
and the known arbitrary initial state ^-s 

\m)^e-^"'\m)- (E.27) 

In coordinate representation g{x, x'; t) = {x\ exp{—iHt)\x') is referred to as the retarded 
Green function. It may be written, using Laplace transform techniques, in terms of the 
outgoing Green function to the problem, 

x'; t) = ^ j G+{x, x'; fc)e-'''='*/2" 2Mfc, (E.28) 

Co 

where Co stands for the Bromwich contour and corresponds to an hyperbolic contour 
along the first quadrant on the complex k plane. 

One may write a representation of the time-dependent Green function involving reso- 
nant states plus an integral contribution in a similar fashion as discussed for the half-line 
in [280] . This may be done by deforming appropriately the integration contour Co on 
the complex fc-plane. Since the variation of G^{x,x';k) with k is at most exponential, 
the behavior of the integrand with k in Eq. (jE.28p is dominated by exp{—ihk'^t/2m). A 
convenient choice is to deform it to a contour involving two semi-circles Cs along the 
second and fourth quadrants of the fc-plane plus a straight line Cl that passes through 
the origin at 45° off the real A:-axis. In doing that one passes over some poles of the 
outgoing Green function. In general these include bound states and the subset of com- 
plex poles associated with the so called proper resonant states, i.e., Hekp > Imkp . By 
extending the integration contour up to infinity one obtains that the semi-circles Cs yield 
a vanishing contribution so that one is left with a infinite sum of terms plus an integral 
contribution [280] . namely, 

oo 

g{x,x';t)^ Y,un{x)un{x')e-'''"'/'' 

n=l 

G+{x, x'; fc)e-»'*'='*/2™ 2kdk. (E.29) 

Cl 

Here the summation includes all bound states Ub, and all proper resonant states, Up. In 
this last case, the energies are complex, i.e., Ep ~ £p — iTp/2. In order to derive the 
above expression, one makes use of Eq. (|E.8p for the residue of a complex pole, with 
the resonant functions normalized according to the condition given by Eq. (jE.9p . One 
may manipulate the integrand in the integral term of Eq. (IE.29P , as in Ref . J80] for the 
half-line, to write the retarded Green function as a sum over bound and resonance terms 
plus a continuum integral term along the negative imaginary energy axis, a form that 
resembles the usual expansion as a sum of bound states and a continuum contribution 
along the real energy axis. 



^ (2m) 27r / 
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The time-dependent expansion given by Eq. (|E.29p holds for arbitrary values of x and 
x' . However, by restricting these values according to the condition ^ (a;, x') ■< L, 
discussed above, one may obtain a time-dependent expansion that includes the full set 
of bound, antibound and resonant states (associated with both the third and fourth 
quadrant poles of the complex /c-plane). Noticing that 

1 _ 1 

allows to write Eq. (jE.15|l . using Eq. (jE.16p . as 

The above representation is very convenient because substituting it into the integral term 
m Eq. ([R29l) leads to an expansion over the full set of bound, antibound and resonance 
states of the system, similar to that obtained for the half-line [287l288j , 

oo 

gix,x';t)= u„(x)u„(x')M(fc„,i), (E.32) 

n— — oo 

where M{kn, t) = A/(0, fc„, t) stands for the Moshinsky function. 



1 



(E.30) 
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